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SUMMARY 


The objective of the investigation was to develop the ana- 
lytical tools (analysis and computer simulation) needed to explain 
and predict the dynamic operation of air cushion landing systems 
(ACLS) . The work was carried out through the five tasks listed 
below. 

1. Development of improved analytical models for the fan 
and trunk 

2. Formulation of a heave-pitch-roll analysis for the 
complete ACLS 

3. Development of a general purpose computer simulation 
to evaluate landing and taxi performance of an ACLS 
aircraft 

4. Verification and refinement of the analysis by compari- 
son with test data obtained through lab testing of a 
prototype cushion 

5. Demonstration of simulation capabilities through typical 
landing and taxi simulations of an ACLS aircraft. 

Initial results show that fan dynamics have a major effect 
on system performance. Comparison with lab test data (zero 
forward speed) indicates that the analysis can predict most of 
the key static and dynamic parameters (pressure, deflection, 
acceleration, etc.) within a margin of 10-25 percent. The 
simulation can thus be used to evaluate existing configurations 
and develop new designs. 

The computer program and the user's manual are available from 
Computer Software Management and Information Center (COSMIC) , 112 
Barrow Hall, University of Georgia, Athens, GA 30602, and are 
identified as LAR-12303, "Air Cushion Landing System Program FMA4 
Static and Dynamic Simulation of Heave-Pitch- Roll Motion." 



INTRODUCTION 


Background 

The Air Cushion Landing System (ACLS) concept has been applied 
to several aircraft, such as the Buffalo (XC-8A) , the LA-4 and the 
Jindivik drone. Tests with these aircraft and small scale models 
have demonstrated both the feasibility of the concept and the ad- 
vantages of using such a system in place of conventional landing 
gear. In several instances, however, current designs have exhi- 
bited unexpected behavior, and the need for a better understanding 
of ACLS operation has thus been felt. Recognizing this need, NASA 
initiated a phased research program to develop the analytical tools 
necessary to understand and predict ACLS performance. Once the 
analytical models and verifying experiments are complete, this 
work will then lead to the design and fabrication of improved sec- 
ond generation air cushion landing systems. 

This report represents the completion of the first stage of 
the research program in which the basic analysis has been devel- 
oped and subjected to initial verification. The analytical model, 
based on the fundamentals of solid body and fluid mechanics, deter- 
mines the heave, pitch and roll motion and pressure, flow and 
acceleration of an ACLS during landing impact, slideout and taxi. 
The model has been verified by comparing the analytical results 
with those obtained from controlled static and dynamic experiments 
with a prototype cushion in the laboratory. A computer program 
which incorporates the ACLS analysis has also been developed to 
evaluate current ACLS designs, and to help in the development of 
improved configurations. 

Subsequent sections of this report describe details of the 
analytical model, the results of the experimental verification, 
and a typical application in which the analysis has been used to 
evaluate the characteristics of the Buffalo aircraft. The computer 
program itself is not described in this report, but is available 
to interested users through NASA COSMIC services (identification 
LAR-12303) . 
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State-of-the-Art Review 


The feasibility of Air Cushion Landing Systems was estab- 
lished when an ACLS LA-4 aircraft first flew in 1967. Since that 
time, several aircraft applications have been studied by the Air 
Force, Navy and NASA (see Table I) and two of these (Buffalo and 
Jindivik) have been developed to the point of full-scale testing. 
Scale model tests of cushions and entire aircraft have also been 
carried out, including static tests, drop tests, forward speed 
tests and wind tunnel tests. 

Over the last ten years, a number of reports and papers have 
been published describing this work, from initial feasibility and 
design studies to subsequent analyses, computer simulations and 
testing. This literature is summarized in Table II. 

The first analysis of a complete ACLS was carried out by 
Earl (ref. 1) . His report on the feasibility of using an ACLS for 
the C-119 aircraft included a static analysis of the cushion and 
its subsystems. The static theory was supported by scale model 
tests . 

Digges (ref. 11) also analyzed the complete ACLS, starting 
with an evaluation of the various jet theories which model the 
outflow from the trunk and cushion. He compared the flow para- 
meters (e.g. , cushion-to-trunk pressure ratio) predicted by the 
theories for different jet heights. The jet analysis was then 
expanded to include distributed flow due to distributed trunk holes 
and the results thus obtained were compared with those found from 
the flow restrictor theory, which models the air gap as an orifice 
between the cushion chamber and atmosphere (Figure la, lb) . While 
test results showed that the flow restrictor theory predicted the 
cushion flow more accurately than the distributed jet theory, it 
required specially measured orifice discharge coefficients. 

Digges also analyzed the trunk shape, assuming that the trunk 
material behaves like an ideal membrane (Figure 2a) . Trunk shapes 
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TABLE I 


ACLS AIRCRAFT DEVELOPMENT SUMMARY 


Aircraft 

Type 

Designation 

Sponsoring 

Agency 

Development Level 

Design 

Study 

Model 

Tests 

Full 

Scale 

Tests 

1. Light Amphibian 

LA- 4 

Bell 

X 

X 

X 

2. Transport 






Light 

Buffalo 

Air Force 

X 

X 

X 

Medium 

C-119 

Air Force 

X 

X 


Heavy 

Span Loader 

NASA 

X 



3. RPV 







Jindivik 

Air Force 

X 

X 

X 


Ryan 147G 

Air Force 

X 



4. Fighter 







A- 4 

ARP A/Navy 

X 

X 



F-8 

ARPA/Navy 

X 

X 


5. Space Shuttle 


NASA 

X 







TABLE II 


SUMMARY OF PRINCIPAL ACLS LITERATURE 


Performing 

Do to 

Principal 
Author 6 

Aircraft 

ACLS 

Subsystem 

Anal 

Type 

ysls 

Mode 


Cxpe 

inent 

Forward { 

Hind 

Rcnerke 

AgiTicy 


Reference 

Selected 

Studied 

Static 

Dynamic 

Heave 

Pitch 

Roll 

static Test 

Drop Teet 

Speed Teet 

Tunnel Test ' 


noil 

5/67 

Eorl 

C>119 

Complete 

X 


X 



1/3 Scale 

1/3 Scale 

1/3 Scale 

1/12 Scale j 

Pcanibllity Study 

AurOKpacu 


(1) 


ACLS 






1/12 Scale 





Co. 

8/Cd 

Eorl 

C-119 

Complete 






1/3 Scale i 

1/3 Scale 


Extension of 

1/3 Scale 



(2) 


ACLS 







1/4 Scale 


Earlier Work 


9/69 

Stauffer 

LA-4 

Complete 










Prototype ACLS 



(3) 


ACLS 





X 





Aircraft 


12/70 

Ryken 

Space 

Complete 

X 

X 

X 

X 









(4) 

Shuttle 

ACLS 












‘j/V2 

Cuius 

liui'falo 

Cotupluto 


* 

X 

X 


2D Trunk 

1/4 Scale 1 

1/4 Scale , 

1/10 Scale 

Includes Computer 



(5) 


ACLS 






1/4 Scale 

1/10 Scale 

l/IO Scale 


Simulation 











1/10 scale 




Dovnlounicnt 


7/72 

Ryken 

Ryan 

Complete 



X 







Design Study 



(6) 

147G RPV 

ACLS 











3/7< 

Rykon 

Jindivik 

Complete 

X 

X 

X 

X 






Similar to Coles* 



(7) 


ACLS 










Analysis 


S/74 

Eorl 

Boeing 

Complete 

X 


X 







Design Study 



(6) 

737 

ACLS 











11/75 

Earl 

LA- 4 

Suction 

X 












(9) 


Braking 











Wright 

1970 

Goman 

. 

Peripheral 

X 





2D Trunk 





tactoraon 


(lOJ 


Jot Airflow 











AFB 

6/71 

Oigges 

_ 

Complete 

X 

X 

X 



2d Trunk 6 

Prototype 



Computer Simulation 



(11) 


ACLS 






Prototype 

Cushion 

Cuehlon 



Development 


0/72 

Vaughan 

Jindivik 

Complete 

X 

X 

X 







Computer Slnulation 



(12) 


ACLS 










Development 


9/72 

Voughon 

Buffalo 

Complete 













(13) 


ACLS 






X/4 Scale 

1/4 Scale 





1973 


Jindivik 







Ejector and 







(14) 

Source 






tip turbine 
fan 






1973 


Jindivik 







1/10 Scale 

1/10 Scale 

1/10 Scale 


Results Compared with 



(IS) 

ACLS 










Full Scale Model 


9/73 

Rogers 


Pcriohcral 

X 





2D Trunk 







(16) 


Jet Airflow 












9/73 

Rodrigues 

Buffalo 

Complete 






1/10 Scale 

1/10 Scale 



Results Compared with 



(17) 


ACLS 










X/4 scale 


3/74 

Carreras 

- 

Trunk/Cushion 

X 





2D Trunk 







(1(0 


Airflow Syatem 












3/74 

Daucr 

Jindivik 

Mechanical 

X 

X 

X 

X 

X 





Lumped parameter 



(19) 


System 










nocnanical model 


3/74 

Stuart 

Buffolo 

Complete 

X 


X 







Model Developed from 



(20) 

(Scale 

Kodcls) 

ACLS 










previous expur^ments 


4/74 

Vaughan 

Jindivik 

Complete 






Full Sca\c 

Full Scale 



Includes pull tests 



(21) 


ACLS 












<n 


TABLE II (CONCLUDED) 

SUMMARY OF PRINCIPAL ACLS LITERATURE 


Performing 

Acjer.cy 

Date 

Principal 
Author k 
Rcforciico 

Aircraft 

Seloctod 

ACLS 

Subsyuten 

Studied 

Static 

Ana 

Type 

Dynamic 

lysis 

Mo 

Heave 

do 

Pitch 

Roll 

Static Test 

Exp 

Drop Test 

:riment 
Forward 
Speed Teat 

Hind 

Tunnel Test 

Remarks 

Ohio State 

5/71 

Han 

(22) 

- 

Cuuhion 

Airflow 

X 


X 









3/72 

7/75 

Han 

(23) 

Han 

(24) 

- 

Cushion 

Airflow 

Cushion 

Airflow 






2d Cushion 
Airflow test 
apparatus 
2D Cushion 
Airflow test 




Experimental 
Corroboration of 
Earlier Theory 


1975 

Forxono 

(25) 

Jindivlk 

ACLS Trunk 






Trunk Flutter 
Test Apparatus 




Trunk Flutter 
Experiments 

Foster- 

Miller 

Associates 

5/75 

Captain 

(26) 

- 

Complete 

ACLS 

X 

X 

X 

X 






Computer Simulation 
Development 

N'SSRC 

2/73 

Lee 

(27) 

A-4 

Aerodynamic 

Configuration 









22t Model 


Boeing 

3/73 

Milne* 

A-4 (1/3 
Scalo) 

Complete 

ACLS 

X 


X 



1/J Scale 

1/3 Seale 




NASA 

9/73 

Thompson 

(29) 

C-8 

Complete 

ACLS 






1/LO Scale 


1/10 Scale 




3/75 

Lelar.d 

(30) 

Navy 

Fighter 

Complete 

ACLS 






0.3 Scale 

0.3 Scale 

0.3 Scale 



B.r. 

Goodrich 

7/7 4 

Ho** 

Jindivik 

Complete 

ACT.S 

X 

X 

X 

X 






Computer Simulation 
OovoLopment 

Soot.vwesL 

Research 

Institute 

6/73 

Bass 

(32) 

•• 

ACLS Trunk 






Trunk 

Flutter Tost 
Ai^paratus 




Trunk Flutter 
Experiments 

Sendair 

1/73 

KeCuddon 

(33) 

Jindivik 

Complete 

ACLS 

X 









Design Study 

Lockheed 

12/72 

Hurkamp 

(28) 

Span 

Loader 

Complete 

ACLS 

X 









Design Study 


MUns, Petrr; and Gardner, Lloyd; A-4 Twin Pod Air Cushion Test 
Specimen. Report D180-15370-1, tiie Boeing Company, Seattle, 
Washington, March 1973 

Ho, P.H., Ajialysis oC the Transient Dynamic Pitcl» Characteristics 
of i»n Mr Cushion Rccovory System ;'or the JindiviK Aircraft. 

B. r. Coodrich Research and Development Center, Ohio, July 19/4 






Figure 1. Analytical Models for Cushion and Trunk Flow 





(b) Membrane Trunk - With Hoop Stiffness (ref. 5 ) 



(c) "Frozen" trunk (ref. 26 ) 

Figure 2. Analytical Models for Trunk Deformation 
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were predicted for various cushion-to-trunk pressure ratios and 
compared with the test results obtained with a straight (2 dimen- 
sional) section of trunk. This membrane trunk analysis however 
did not include the effects of the hoop stress and deflection of 
an actual (3 dimensional) toroidal trunk.* The trunk shape and 
cushion flow analysis was used to develop a dynamic model of the 
ACLS. A dynamic simulation computer program was also developed. 
Initial test data from a prototype cushion (at zero cushion pres- 
sure) were found to be in general agreement with the computer pre- 
dictions . 

A vertical energy absorption analysis for ACLS was performed 
by Vaughan (ref. 12) , who also developed a dynamic simulation com- 
puter program. This program was restricted to heave mode simula- 
tion. A parametric analysis was performed to determine the effects 
of the various parameters on the load-stroke response of the sys- 
tem. 

An analysis of pitch motion of the ACLS is described by Coles 
(ref. 5) . In this analysis the trunk is divided into sections and 
the forces and moments for coupled heave-pitch motion are obtained 
by analyzing each section separately and then summing the results 
for the whole cushion. Hoop tension effects, not considered by 
Digges, were approximated in this study by introducing an experi- 
mentally calculated radial stiffness in the force balance analysis 
of the trunk by means of an external spring (Figure 2b) . Trunk 
dynamics were modeled by using auxiliary pressure variables re- 
lated to the actual pressure variables by first order lags. A 
computer simulation, which included aircraft aerodynamics, was 
developed to generate heave-pitch motion predictions of the ACLS 
during landing impact. Ryken (ref. 7) expanded Coles' analysis 


* 

In this report, the term hoop refers to the trunk lines that run 
around the (oval) perimeter of the cushion and trunk. 
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by introducing two additional system pressures (the forward gap 
and the aft gap) calculated from an experimentally observed expo- 
nential relationship. The analysis was used to update the com- 
puter program developed by Coles, and various simulations were 
carried out to predict the behavior of the Jindivik. 

In an attempt to include analytically the effects of circum- 
ferential tension. Captain (ref. 26) assumed the trunk to be 
"frozen", i.e., the trunk shape outside the ground contact zone 
to be fixed, and independent of pressure (Figure 2c) . With this 
assumption, a dynamic heave analysis was developed and subsequent- 
ly extended to include pitch motion. Computer simulations were 
developed for both the heave and heave-pitch modes. 

A dynamic heave-pitch simulation of the Jindivik is described 
by Ho* who solved existing analytical models (refs. 2, 23) on a 
hand calculator. 

Simpler approaches to modeling the ACLS have also been taken. 
In a preliminary study by Bauer (ref. 19) , the ACLS is represented 
by a spring-mass-dashpot model, and the various system parameters 
are quantified from test data obtained with the full scale 
Jindivik. Stuart (ref. 20) proposed a similar type of model in 
which several of the equations were quantified through experiments 
with scale model cushions for the Buffalo. 

In addition to studying the entire ACLS, several reports deal- 
ing with specific ACLS subsystems have also been published. The 
majority of these reports have dealt with the modeling of the fluid 
interaction between the trunk, cushion and atmosphere. Most of 
these studies have included experiments on two dimensional trunks 
with various orifice configurations. Theories used to predict the 
experimental behavior include the jet theories such as the Barratt, 
thin jet and exponential theories (Figure la) and the flow restric- 
tor theory (Figure lb) . The jet theories generally use change of 
jet momentum to predict the cushion pressure, whereas the flow re- 
strictor theory considers the gap area under the trunk as an 


See Table II for reference 
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orifice and predicts the cushion pressure from the orifice equation 
(refs. 11, 16) . The general conclusion drawn from these studies is 
that the orifice theory is more suitable for ACLS than the jet 
theories, apparently because a discrete jet pattern is not estab- 
lished. Han (ref. 22) has obtained theoretical results by approxi- 
mating the trunk-cushion interaction by the airflow through a 
porous plate (Figure Ic) . An experimental verification is also 
described (ref. 23) . The theory works well for low gap heights 
characteristic of trunk-ground contact, but is not attractive for 
the larger clearances. More recently, Han (ref. 24) has reported 
some experimental work on a distributed jet emanating from single, 
double and triple slot configurations. About the only published 
studies on ACLS fans is the work of Kunstadt (ref. 14) who evalua- 
ted both ejectors and tip turbine fans for the Jindivik. 

The general conclusions that emerge from the state-of-the- 
art review of ACLS model development are as follows. 

1. Operation of the trunk-cushion-atmosphere airflow system 
should be modeled by an orifice theory rather than a jet 
theory. 

2. The trunk model is a known bottleneck that has restric- 
ted the advancement of ACLS analysis. Of the three 
trunk models presently available, the membrane model 
(Digges) is too "soft” because it neglects hoop tension 
effects, the "frozen" trunk model (Captain) is too 
stiff, because it overcompensates for hoop tension, and 
the model developed by Coles is confined to modes and 
applications where special trunk test data are available 
to quantify the model. 

3 . The pressure distribution under the trunk in the ground 
contact zone has not been conclusively established. 

Most analyses assume that during ground contact the 
trunk forms a perfect seal with contact pressure equal 
to the trunk pressure. Very little data are available 
to validate this assumption for a real three dimensional 
trunk, and further analysis and testing is needed to 
establish the model for trunk contact pressure and 
sealing . 
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4. 


The air source has not received much attention, and 
the effects of source dynamics on ACLS landing behavior 
are largely unknown. 

5 . Although several computer models have been described in 
the literature, there is no generally available and 
experimentally verified simulation that can be used 
without special test data to evaluate and improve ACLS 
designs. 

The conclusions of the state-of-the-art review helped define 
the specific tasks needed to address those critical areas that are 
as yet unresolved. The work was broken down into five tasks: 


1. Development of improved models for the fan and trunk 

2. Formulation of a heave-pitch-roll analysis for the 
complete ACLS 

3. Development of a general purpose computer simulation 
to evaluate landing and taxi performance of an ACLS 
aircraft 

4. Verification and refinement of the analysis by compari- 
son with test data obtained through lab testing of a 
prototype cushion 

5. Demonstration of simulation capabilities through typi- 
cal landing and taxi simulations of an ACLS aircraft 

The preceding section has summarized previous work in the 
field and identified the key areas for further study. The static 
and dynamic analysis of the ACLS is developed in the next section. 
All major assumptions are discussed, and the complete ACLS analyti- 
cal model, including improved models for the trunk and fan, is 
derived. The full set of equations that describe the model are 
presented in the Appendix. The test setup and the static and 
dynamic experiments performed on the prototype cushion* are then 
presented. The experiments carried out include static heave, 
pitch and roll tests, in which various ACLS parameters, including 

*One half of a 1/3 scale dual cushion landing system built by 
Boeing under an earlier Navy contract. 
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ground contact pressure distribution, were measured for different 
loading conditions. The dynamic tests consist of drop tests from 
different heights with different angles of attack, pitch and roll 
moment release tests, and a coupled heave-pitch start-up test. 

The validity of the analytical model is evaluated by comparison 
with the test results. Finally, the capabilities of the analyti- 
cal model are demonstrated by carrying out a computer simulation 
of the Buffalo aircraft. The results presented include heave and 
pitch motions, pressures and accelerations during a typical 41 m/s 
(80 knot) approach and landing, and during taxi over an ungraded 
runway . 

Names of several manufacturers and their products have been 
referred to in the report. This does not constitute official 
endorsement of such products or manufacturers, either expressed or 
implied, by the National Aeronautics and Space Administration. 

ANALYSIS 

Basic Configuration 

The basic ACLS configuration analyzed is shown in Figure 3. 

The model includes four primary subsystems: the fan, the feeding 

system, the trunk and the cushion. The configuration of these 
systems has been chosen sufficiently general so that they can 
represent a variety of practical designs. Air from the fan flows 
through the ducts and plenum (feeding system) and enters the trunk. 
The trunk has several rows of orifices that communicate both with 
the cushion and. atmosphere. Thus, the airflow from the trunk has 
two components - one part entering the cushion and the other 
leading directly to the atmosphere. The cushion flow exhausts to 
the atmosphere through the clearance gap formed between the trunk 
and ground. In addition to the basic flows described above, two 
other flows have been included in the model for generality. These 
are the plenum bleed flow and the direct cushion flow. Plenum 
bleeding causes some of the air to flow directly from plenum to 
atmosphere, and has been used in some designs to improve the dy- 
namic characteristics of the air supply system. Direct flow from 
the plenum to the cushion can also improve dynamic response. A 
pressure relief valve is also included in the basic configuration. 
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II 


It allows additional flow to vent from the plenum whenever the 
pressure exceeds a preset level, and thus improves stability by 
reducing fan stall. 

The support force acting on the aircraft is made up of two 
components. The first occurs due to the cushion pressure acting 
over the cushion area. The second, which comes about only during 
ground contact, is given by the contact pressure acting over the 
trunk contact area. The support force, in general, also gives 
rise to a moment, given by the product of the force and its dis- 
stance from the CG of the aircraft. 

In plan, the cushion has an oval shape, made up of a rectan- 
gular section with semicircular ends. The lengths a and b are 
the horizontal and vertical spacing between the points of attach- 
ment of the trunk to the aircraft body. The initial (undeformed) 
trunk shape is defined by the above two parameters, and the peri- 
meter i and height Hy as shown. is the (uniform) spacing 

between the rows of peripherally distributed orifices. The number 
and orientation of the orifices is selected independently by the 
number of orifice rows N^, the number of orifices per row and 

the orientation parameter a . The cushion volume consists of 

P 

two parts: an active (dynamically varying) region and a dead 
(static) region. The active volume depends on the trunk shape 
and ground profile. The dead volume, which is a design variable, 
includes recesses in the cushion cavity as shown. 
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Method of Approach 


The configuration of the cushion was described in the last 
section. Here, a description of the overall system and the initial 
steps taken to set up the analytical model are presented. 

The landing system forces transferred to the aircraft act 
through the cushion and trunk. To help calculate these forces, the 
trunk and cushion are divided into segments as shown in Figure 4 . 
Each straight section of the cushion and trunk is divided into 2M 
rectangular segments, while each curved end is divided into 2N pie- 
shaped segments. Thus the total number of segments is 4 (M + N) . 

The cushion and trunk parameters are calculated first for each seg- 
ment and then summed to give their total system values. 

The dynamic analysis of the vehicle system is best derived 
with the help of two orthogonal coordinate frames of reference: a 
coordinate frame fixed in space (inertial frame) , and a coordinate 
frame fixed to the vehicle (vehicle frame) with origin at the 
aircraft CG. The reason for two frames can be appreciated by 
recognizing that 

(a) Newton's law for translational motion can be 
directly applied to the CG acceleration expressed 
relative to the inertial frame. 

(b) The corresponding law for rotational motion can be 
applied more conveniently in the vehicle frame, be- 
cause rotational inertia about any vehicle axis is 
constant, while the rotational inertia about any 
inertial (fixed) axis varies with aircraft position. 

Accordingly, the two frames of reference have been defined 
as shown in Figure 5. The vehicle frame with origin at the air- 
craft CG has roll, yaw and pitch axes x, y and z respectively 
fixed to the aircraft body as shown. The inertial frame has cor- 
responding axes X, Y and Z fixed in space. The two frames coin- 
cide only when the aircraft has not undergone any rotation from 
equilibrium. 
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Figure 4. Division of Trunk into Segments 
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In the analysis, the actual runway profile underneath the 
ACLS is approximated by segments that coincide in plan with those 
of the trunk and are parallel to the cushion hard surface as shown 
in Figure 5. With this model, all pressure forces act parallel 
to the vehicle yaw axis so that the segment torque component 
about the aircraft CG can be easily computed by multiplying the 
segment force by the fore-and-aft and/or lateral separation be- 
tween the segment and the CG. 

The model sketched in Figure 5, and the analysis derived sub- 
sequently, has been set up to cover the heave, pitch and roll mo- 
tion of the aircraft during touchdown, slideout and taxi over an 
irregular runway. Sideslip forces and yaw moments about the 
vehicle axis have not been taken into account at this time on the 
ground that, in practice, these will be determined by the response 
of the as yet undefined aircraft steering system. The forward 
velocity of the aircraft, determined independently of the dynamics, 
is found from the initial aircraft velocity and the selected 
braking coefficient. 

As mentioned earlier, aircraft translational motion is found 
by applying Newton's law to the aircraft CG in the inertial frame. 
Thus for heave (vertical motion) 


Y = F 
a eg y 


( 1 ) 


where 


M is the aircraft mass 
a 


Y is the aircraft CG motion in the inertial 
vertical direction 


and 



is the ACLS force in the inertial vertical 
direction 
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Similarly, for forward motion 


where 


X = -u 
eg 


X is the aircraft CG motion in the inertial 
forward direction 


(2) 


and 

y 

is the 

For 

aircraft 

rotation 

gives 

-y 

T 

A 

dt 

where 

T 

is the 

and 

H 

is the 


is the braking coefficient. 


H 

is the torque vector about the aircraft CG 
is the angular momentum vector. 


(3) 


Substituting for the generalized angular momentum vector and 
expressing the above equation in matrix form gives 


1 

1-3 

X 

1 
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• CD 

1 




1 

• CD 

1 

T 

y 

= 

I 

• 

ip 

+ 


X 

I 

•• 


1 

N 

1 




- <P - 


4 . 


_ _ 


. (p . 


where T , T and T are the torque components acting about the 
X y z 

vehicle x, y and z axes 

• • •• •• 

0, ip , 9, 4^, and (J) are the roll, yaw and pitch velocities 

and accelerations along the vehicle x, y and z axes 
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and the inertia matrix [l] is given by 




r I 

I 

I 1 



X 

xy 

zx 

I 

= 

I 

I 

I 



xy 

y 

yz 



I 

I 

I 

La <J 


L zx 

yz 

z J 


where I , I , I are the moments of inertia of the aircraft along 
X y z 

the vehicle x, y and z axes. 


and 



I are the cross products of inertia. 

Z X 


The translational equations of motion, Eqns (1) and (2) , can 
be easily integrated to give the X and Y coordinates of the 
aircraft CG as a function of time. Integration of the rotational 
equation of motions, Eqn (4) , however is less straightforward, be- 
cause the end result of the integration which is angular rotation 
is not a vector and thus cannot be uniquely defined by the angles 
0, and <j) as obtained from the solution of Eqn (4). In problems 
of this sort, in which it is necessary to uniquely determine the 
orientation of a body rotating about an arbitrary (and moving) axis, 
the difficulty is avoided by a coordinate transformation which 
expresses the rotational equations in terms of Euler angles. 

Euler angles signify nothing more than angular rotations about 
prescribed axes oarv-ted out in a defined sequence . Unlike conven-- 
tional angles, the definition of a specific sequence for adding 
Euler angles allows them to uniquely fix the orientation of a body 
undergoing rotation about an arbitrary axis. 
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The Euler angles and rotation sequence adopted for* the subse- 
quent analysis is shown in Figure 6. The three steps needed to 
fix the angular orientation of the aircraft are as follows: 

1. Beginning with the vehicle axes parallel to the inertial 

axes, the first Euler angle rotation is made about 

the inertial yaw axis Y which, at this stage, coin- 
cides with the vehicle yaw axis, y. Note that after 
making this rotation, the inertial and vehicle pitch 
and roll axes are no longer parallel. 

2. The second Euler angle rotation 0 is made about the 
vehicle roll axis x. Note that after this rotation, 
the inertial and vehicle yaw axes are no longer parallel, 
and the orientation between vehicle and inertial pitch 
axes has further changed. 

3. The third Euler angle rotation is made about the 

vehicle pitch axis z. Note that after this rotation, 
the vehicle roll axis has been displaced from the posi- 
tion it had during Step 2. 

Thus, in summary, the first rotation is about the inertial 
yaw axis and the third rotation is about the final pitch axis; 
both of which are easy to identify. The second rotation however 
is about the -intermediate vehicle roll axis, which is different 
both from the inertial roll axis due to the rotation of Step 1, 
and from the final vehicle roll axis due to the rotation of 
Step 3. Of course, when only pitch-heave simulations are carried 
out (i.e., when roll excitation is not present), then the Euler 
roll and yaw angles are zero = 0^ = O) and the Euler pitch 

angle cj)^ coincides with the usual inertial pitch angle ip. 

With the above definition, the coordinate transformation can 
be carried out to give the relationship between the angular velo- 
cities in the vehicle frame (0, ifi, 4>) and the corresponding Euler 
angles and their derivatives. The derivation, which is presented 
in the Appendix shows that: 
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0 


(5) 


= i sin d) cos 6+0 cos cb 
e e e e 


ill = il) cos d) cos 0 — 6 sin A (6) 


^ sin 6^ (7) 

Relationships for vehicle frame acceleration (0, rp and can 
also be obtained by differentiating the above equations with respect 
to time. When these relationships are substituted into Eqn (4) , 
the differential equation can then be solved to give the Euler 
angle components, and thus fix the angular orientation of the air- 
craft as a function of time. 


Analytical Development 


The analytical model of the ACLS consists of a set of equa- 
tions which when solved determines the pressures, flows, forces 
and motion of the system as a function of time for various air- 
craft and runway parameters. There are two parts to this model: 


(a) The Static Model s which comprises the relationships that 
determine the static characteristics of the system at 
equilibrium. In addition to providing design data these 
relationships also provide the initial conditions for 
the dynamic model. 

(b) The Dynamio Model s which comprises the differential 
equations of flow and motion (state equations) from 
which the pressures, flows, forces and motion can be 
determined as functions of time. 

Before proceeding with the analytical development, it is 
helpful to divide the overall ACLS into two interrelated systems: 
the flow system and the force system. These systems are shown in 
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Figures 7 and 8. The flow system establishes the pressure-flow 
relationship for the various subsystems of the ACLS. The force 
system establishes the corresponding force-motion relationships . 

The interdependence of the two systems comes about because the 
trunk deflection obtained from the force system changes the vol- 
umes and orifice areas that form part of the flow system. Similar- 
ly, the cushion and trunk pressures found from the flow system 
give rise to forces and moments that form inputs to the force sys- 
tem. With the definition of these two systems, the analytical 
development can now proceed. The primary equations of the static 
and dynamic models are derived in the subsequent sections. The 
complete set of equations and symbol list is given in the Appendix. 

The Static Model 

The equilibrium conditions are found as follows; 

(a) By satisfying the fan flow constraint, 

°fx = * (^f) 

where the functional relationship is found from the 
static fan characteristic. 

(b) By applying the steady-state flow continuity equations 
to the plenum, trunk and cushion cavities (Figure 7) . 


°fx 

“pa 

+ 

Qpt Qpc / Qv 

(9) 

O 

rt 

II 

®tc 

+ 

Qta 

(10) 

Qca = 

pc 

+ 

Qtc 

(11) 
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(c) By satisfying the trunk shape constraint for -the 
trunk sides. 



where the functional relationship is found from the 
trunk model described later. 


(d) From static force balance in the vertical direction 
(Figure 5) 

= (^c ^c ^t ^cn) 


where 



- aircraft mass 

- cushion pressure 

- cushion area 
trunk pressure 

trunk area in ground contact 
Eulerian pitch angle 

- Eulerian roll angle 
acceleration of gravity 


(e) Prom static torque balance in pitch (Figure 5) . 


( 




T cos 4>^ cos e 
nz ^e e 


(14) 


where 


the total cushion and contact pressure torque 
obtained by adding the torque contributed by each 
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segment of the finite element trunk and cushion 
model, and is given by 


nz 


4 (M+N) 

E 

i = 1 






( 15 ) 


where 


A^(i) 






x^^(i) 


4 (M+N) 


- cushion pressure (gage) 

- trunk pressure (gage) 

cushion area of the i^^ segment 

trunk contact area of the i^^ segment 

X axis distance between the centroid 
of the i cushion segment and the 
geometric center of the cushion 

X axis distance between the centroid 
of the i^ trunk contact segment and 
the geometric center of the cushion 

- number of segments in cushion and trunk 
model 
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(f) From static torque balance in roll (Figure 5) 


/m = -T cos 4> cos 0 (16) 

\ J f nx e e 

where is the total cushion and contact pressure torque 

about the vehicle x axis, and is given by 


4 (M+N) ^ 

V = E [''c 

i=l 




z axis distance between the centroid of 
th 

the i cushion segment and the geo- 
metric center of the cushion 

z axis distance between the centroid 
of the i^ trunk contact segment and 
the geometric center of the cushion 

and the other symbols are defined in (e) above. 

The Dynamic Model 

The state equations for the system are derived as follows: 


where ^ch^^^ 
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(a) By satisfying the dynamic fan flow constraint (see 
"The Fan Model") 


dQ 


fx = 


dt 


P. + P ^ - P 

f af £ 


(18) 


where 


is the static fan pressure rise for the flow 
®fx 


af 


is the fan inlet pressure 


Pp is the plenum pressure 


and 


is the fluid inertance of the fan ducts 


(b) By applying the flow continuity equation to the plenum 
trunk and cushion cavities (Figure 7) . The dynamic 
continuity equation requires that the net inflow into 
the cavity equals the rate of increase of fluid mass 
within the cavity. 


^ °pa-V-°PC-«v)o 

^ ('’''t) = (V °to - °ta)<= 

A (V * °tc - °ca) '21) 
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(c) By satisfying the trunk shape constraint for the 
trunk sides, 

(Vv) “ 

(d) From the dynamic force balance in the vertical direction 
(Figure 5) 

MY =/f +F. + Fj^ + F . \cos (j) cos 9 

a eg \ cp tp df ct) e 

- M^g (22) 


where 

o 

t! 

1 

cushion pressure force 


tp 

trunk contact pressure force 


F - 

df 

aerodynamic drag force 


o 

rt 

1 

trunk damping force 

These 

forces are 

evaluated as follows: 


F 

cp 

P A 
c c 


F. 

tp 

^t ^cn 


F = 

^df 

" ^ph P 


where V is the heave velocity component. (Trunk damping 

force is evaluated in the trunk model described 
later. ) 
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where 


(e) By applying the moment of momentum theorem about 
the vehicle pitch axis (Eqn (4) , with 
^ = ip = 0 and T = o) 


$ = i fT - e(i e + i ] 

X I z \ xy yz ^ / J/ \ X z xz j 

-I Ft + 4) /le + i I - \ 

zx [ X ^ \ xy yz jJA X z xz ; 


( 23 ) 


T=T '+T +T +1^ 

X ■'cpx -^tpx ^dfx ^tx 


(24) 


cpx 


’tpx 


- total torque on vehicle x-axis 

- cushion pressure torque on vehicle 
x-axis 

- contact torque on vehicle x-axis 


dfx 


torque due to aerodynamic force on 
x-axis 


tx 


- torque due to trunk damping force on 
vehicle x-axis 


(f) By applying the moment of momentum theorem about 
the vehicle roll axis ^Egn (4) with 
ip = ip = 0 and ~ ^ | 

0 = [t + 4 /i e + I I - I 

z [ X \ xy yz ^/J/V z x zx / 


- ".X [’’z - ® Ky ® <25) 
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where 


where 


is the total torque on the vehicle z axis whose 
components are defined analogous to those of T in 
(e) above. 

(g) By carrying out a coordinate transformation and 

expressing the vehicle frame velocity and accelera- 
tion components in terms of Euler angles and their 
derivatives as discussed earlier. 


0 = 

sin (j) 

e 

cos 

0 + 
e 

0 cos 
e 

-e- 

(D 

(5) 


cos d) 

e 

cos 

0 - 

e 

• 

0 sin 
e 

tf) 

e 

(6) 

^ 

• 

- dj sin 0 
^e e 




(7) 


(h) By applying Newton's Law to the pressure relief valve 
which is modeled as a second-order mass-spring- 
dashpot system with preload and nonlinear stops. 

The valve displacement is given by: 

m X + Z X + k X = (p - P,\a (26) 

vv vv vv yp pb / V 

subject to the constraint that limits valve motion 

within the range 0 < x < x . 

^ ~ V “ a 


m^ - relief valve mass 

Zy - valve damping constant 

- valve spring stiffness 
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valve motion 


X 


V 



plenum pressure 

valve preload pessure (cracking pressure) 
valve area 


The valve outlet diameter and discharge coefficient form 
external inputs to the model. 


Subsystem Models 


The Fan Model 

Almost all the work in the past has treated the fan as a static 
element, i.e., as a flow source whose output depends only on the 
instantaneous value of the backpressure. While this model is ade- 
quate for slowly changing pressures, it leads to unrealistically 
fast flow changes when the ACLS experiences touchdown. Clearly, 
during touchdown and other rapidly changing conditions, dynamic 
phenomena that introduce lags in the fan flow response become im- 
portant, and must therefore be included in the fan model. 

There are several levels of fan models that can be developed 
to characterize dynamic fan flow. In this work, a lumped parameter 
model has been formulated, which characterizes the fan as a static 
pressure rise element in series with an inertance and capacitance 
as shown in Figure 9a. Flow from the atmosphere goes through an 
inlet orifice and then experiences a pressure rise in the static 
element as set by the steady-state fan characteristics. In dynamic 
operation, the flow must also experience a pressure change due to 
the air mass that has to be accelerated. This is modeled by the 
stream tube inertance. Finally the compressibility of the air 
changes the air mass contained in the fan passages, and thus affects 
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(a) Dynamic Model 



(b) Test Setup 


Figure 9. Fan Model Development 



the output flow. This characteristic is modeled by the volume 
(capacitance) of the fan passages. Thus, during static operation 
or when the pressure is changing very slowly, the output flow 
equals that found from the static fan curve. However, when the 
pressure changes suddenly, the air mass and compressibility intro- 
duce lags in the system and prevent instantaneous flow changes 
from taking place. 

Two simplifications that can be carried out with this model 
are to add the fan volume (capacitance) to the plenum volume, and 
to increse the stream tube inertance to take into account any long 
ducts connecting the fan to the cushion. With these changes, the 
fan flow into the plenum is given by 


■f V dt / 


P ^ - P 

af f p 


( 18 ) 


where 


is the fan flow 


is the fan inlet pressure 


and 


P^ is the static fan pressure rise 

Pp is the plenum pressure 

is the inertance of the fan and connecting 
ducts . 


The drop across the upstream orifice (if present) is given by 
the conventional orifice flow equation 


'“af = P 


(Ofx/Vdu)^/ 


( 27 ) 
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where A is the orifice area 

u 

p is the air density 
and is the discharge coefficient 

The static pressure rise is found from the fan curve, which 
can be represented analytically by a polynomial fit of the form 

Ofx = “O “l “2 + “3 “4 <28) 

The inertance is given by 



where and 5-^^ are the lengths of the fan flow pass- 

ages and the fan-to-plenum duct respectively, and 
and A^^ are their corresponding cross-sectional 
areas . 

The above model is set up so that it can represent a variety 
of fans by using the appropriate characteristic for the static 
pressure rise. No other test data or empirical relationships are 
needed to quantify the model, because the inertance (and capacitance) 
are determined exclusively from fan geometry. 

In parallel with the analytical development, dynamic fan tests 

3 

were carried out by NASA on a 1 m /s (2000 cfm) 5 kPa (100 psf) 

Joy axial flow fan (ref. 31). This fan was chosen because it powers 
the prototype cushion subsequently tested to verify the complete 

i 

model. The test setup is shown in Figure 9b. It consists of the fan 
attached to a duct with a rotating damper (variable orifice) at the 
other end. Pressure taps and an orifice are included to measure 
static flow characteristics, and pressure taps on the fan were 
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used to measure the dynamic pressure-flow characteristic as the 
damper rotates at various speeds. Figure 10 shows the test data 
obtained when the damper opening was changed in increments . The 
high scatter in the left half of the curve occurs because flow 
stability in this region is low, and pressure surges do not readily 
damp out. The test characteristic however is essentially the same 
as that obtained from the manufacturer. The figure also shows the 
polynomial fit\^(Eqn (28)) used to characterize the static element of 
the model. The negative flow characteristic was estimated by 
determining the cross-sectional area of the flow passages and 
assuming that during pressure reversal the fan passages behave as 
orifices with a linear pressure drop characteristic. Figure 11a 
shows the results as the damper is rotated continuously at various 
speeds. The solid line represents damper closing and the broken 
line damper opening. Perturbation frequencies (which correspond 
to twice the damper speed) of 1 Hz and 5 Hz were used because 
they lie within the normal range of ACLS operation. The curves 
show a departure from the static characteristics, particularly at 
the higher frequency. The fan outlet pressure now depends both on 
the magnitude of the flow and its sense (i.e., increasing or de- 
creasing) . Increasing flows result in a lower pressure rise than 
decreasing flows because in the former case, part of the static 
pressure rise is taken up in accelerating the air mass in the fan 
passages, while in the latter case, the air mass is being decelera- 
ted, and so gives up its kinetic energy to increase the pressure 
rise. Previous fan models have been unable to account for this 
behavior. The model developed here overcomes this drawback. Simu- 
lation results using actual values of inertance and capacitance 
(as found from the fan and test fixture drawings) are shown in 
Figure 11b. 

As can be seen, the model predicts the right trends, and 
should therefore result in an improved simulation of touchdown and 
other rapidly changing situations. 
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The Trunk Model 


The trunk model determines the trunk shape parameters (volume, 
and orifice and contact areas ) , contact pressure distribution and 
damping that form inputs to the ACLS flow and force systems. 

Trunk Shape . - In past work, two analytical models have been 
developed for the trunk shape: The Membrane Trunk Model (ref. 11) 

and the Frozen Trunk Model (ref. 26) . (The model of Coles is 
essentially a membrane model.) The shortcoming of both these 
analyses was that they modeled the side and end segments of the 
trunk in the same way while test data now confirm that the shorter 
curved end segments (front and rear) behave very differently from 
the longer, straight side segments. Figure 12 shows the trunk 
cross section measured at the center of the side and end segments 
as the load on the cushion is increased. The entire side segment 
tends to bow outward and avoid ground contact, while the end seg- 
ment remains virtually fixed, except for a flattening in the re- 
gion that actually touches the ground. This difference in behavior 
occurs because the front segment is smaller than the side segment 
and is curved. When the cushion pressure increases due to an 
increase in the load, the radially outward force causes the oval 
trunk planform to become more circular, as shown in Figure 13. 

This causes a hoop tension force T to act around the trunk peri- 
phery. In the side segments, this force acts substantially normal 
to the side excursion 5^ so that its component resisting the 
motion is negligible and the side segment can thus bow outwards 
relatively unrestrained. In the end segments the situation is 
different, since the curvature of the segment causes the hoop ten- 
sion to have a much higher component opposing the motion, so that 
outward motion of the trunk ends is very much smaller t 

Since hoop tension has very little effect on side trunk mo- 
tion, the side segments can be considered as simple two-dimension- 
al membranes, as done in the Membrane Trunk Model. On the other 
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Figure 12. Measured Trunk Profile 
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hand, the fact that hoop tension restrains ("freezes") the trunk 
ends suggests that these segments be modeled by the Frozen Trunk 
Model. Thus the logical step in trunk model improvement is to 
combine the two existing models and form the Hybrid Trunk Model, 
in which the sides are represented by the Membrane Model and the 
ends by the Frozen Model. 

In static heave operation, as the load on the cushion in- 
creases, the trunk ends contact the ground while the sides bulge 
outwards (as shown in Figure 12) to maintain the cushion air gap 
that allows the inflow from the trunk to exhaust to the atmosphere. 
The amount of side trunk bulge (and hence cushion-to-atmosphere 
air gap) adjusts itself such that the resulting cushion pressure 
and cushion-to-trunk pressure ratio satisfies the conditions for 
free (non-contacting) trunk membrane equilibrium and cushion flow 
continuity. In dynamic operation (or in static pitch and roll) a 
situation can exist when the ground clearance is too small for the 
free trunk membrane equilibrium conditions to be satisfied. In 
this situation, the trunk must touch the ground and flatten. When 
modeling this case, only the trunk in the contact zone is consid- 
ered to deform, since trunk-ground friction will tend to prevent 
lateral trunk movement. 

The Hybrid Trunk Model is essentially a limiting case analysis 
of trunk deflection. In general, best results will be obtained at 
the middle of the respective segments, i.e., at the center of the 
side segments, where the trunk behaves very much like an ideal mem- 
brane, and at the center of the end segments, where the trunk shape 
is truly fixed. In the transition region (at and near where the 
segments meet) the trunk will behave somewhere between the membrane 
and frozen trunk approximations. Although an additional refinement 
to model the transition region could be developed, it is appropri- 
ate first to check the results obtained with the basic hybrid 
model, to see whether transition effects are in fact important. 

With this view in mind the overall cushion geometry (height H^ 
and width H^) as predicted by the Hybrid Trunk Model has been 
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compared with test data obtained with the prototype cushion. The 
results are shown in Figure 14 as a function of the Cushion-to- 
trunk pressure ratio pressure ratio increases, 

corresponding to an increased load, the cushion height (which is 
practically equal to the trunk height) reduces, because the side 
trunk lobes move outward, and the end trunk is flattened. When 
the side lobes move out, the cushion width, which is the distance 
between the lowest points on the side trunk segments, increases. 
Because the end trunk is restrained from moving outwards, the 
cushion length remains virtually unchanged. Since the Hybrid 
Trunk Model predicts the actual cushion height and width within 
10 percent, there is little justification at this point for includ- 
ing transition region effects in the model. 

Contact Pressure. - In addition to trunk and cushion shape 
the trunk model also determines the pressure distribution in the 
ground contact zone. The analysis for pressure distribution is 
complicated by the fact that two separate effects must be consid- 
ered; direct trunk-ground contact caused by the trunk pressure 
forcing the trunk against the ground, and airflow through the 
trunk holes into the interstices that remain in the contact zone. 
These two effects are first discussed separately and then combined 
into a single model for trunk contact pressure. 

When two bodies in contact are acted upon by a force F, the 
actual contact occurs at a number of discrete regions rather than 
over the whole area, due to the inherent roughness of the contact- 
ing surfaces. This is illustrated in Figure 15a where the actual 

contact regions (area A^, A 2 A^) are spread over the nominal 

contact area A. The actual pressure distribution (Figure 15b) 

shows high pressure peaks in the contact regions 

and no pressure in the voids. Because the number of contact re- 
gions is large, it is conventional for purposes of a force balance 
to define an average contact pressure P^^ = F/A acting as though 
the bodies were touching uniformly over the entire area A, as shown 
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(c) Average Pressure Profile 


Figure 15. Pressure Distribution between Bodies in Contact 








in Figure 15c. In the event however that one of the bodies is 
porous, flow seepage can still take place through the voids in 
the contact zone, although of course as the average contact pres- 
sure increases, the void region and hence the flow will decrease. 

To help model this flow phenomenon, consider the flow through a 
porous plate held close to, but not touching the ground, as shown 
in Figure 16a. By symmetry, there is a stagnation point at the 
center in which the pressure rises to its original value P^. 

Moving away from the center, the pressure decreases continuously 
until it reaches the ambient pressure P at the edges of the 
plate, as shown in Figure 16b. Thus the driving pressure for flow 
through the ^late at any point is given by the difference between 
the upstream pressure P^ and the pressure in the gap at that 
point. 

For a trunk in ground contact, both discrete contact regions 
and flow through the voids will exist, and a combined force/flow 
model as shown in Figure 17 will be required. This model combines 
the features of the discrete contact and porous plate models de- 
scribed above. The pressure distribution under the trunk in the 
ground contact zone is assiamed to be made up of two components; 
discrete pressure peaks where actual trunk-ground contact takes 
place, and a continuous pressure profile caused by trunk outflow 
through the voids. The assumed profile is shown in Figure 17b. 

For purposes of the force balance, a mean contact pressure acting 
over the nominal contact area can be defined as shown. This mean 
pressure, which is actually the integrated value of the pressure 
profile (including peaks) , can be found very simply by considering 
the equilibrium of a flat segment of trunk in nominal ground con- 
tact (Figure 18) . Isolating a short segment of trunk (area dA) , 
the mean ground pressure is found to equal the trunk pressure, 
since any force T in the trunk itself cannot have a component 
in the vertical direction. Of course, at the very edge of the con- 
tact zone, the change in trunk angle will introduce other terms in 
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(b) Free Body Equilibrium 


Figure 18. Trunk Element Force Balance 




the force balance, but, since this is a localized effect, its in- 
fluence on the mean pressure distribution is negligible. Thus, in 
Figure 17b the mean contact pressure is shown equal to the trunk 
pressure, so that the contact force on the trunk is given by the 
product of the trunk pressure and the nominal contact area. 

For purposes of trunk outflow calculation, the pressure pro- 
file in the non-contacting regions is approximated by a linearly 
decreasing relationship as shown in Figure 17b. The driving pres- 
sure for flow through any trunk hole is thus given by the differ- 
ence between the trunk pressure and the gap pressure at that loca- 
tion. 

As a check of the above hypothesis, tests were carried out to 
measure the base pressure distribution and reaction force of the 
trunk in ground contact. The observed pressure profile, measured 
with manometers at various points in the contact zone, is shown in 
Figure 17c for different loads. Since the manometers only measure 
the fluid pressure and not the local contact presure, only the out- 
flow pressure component is observed. As predicted by the theory, 
this component rises from the cushion pressure at the inner edge, 
to a maximum at the center of the contact zone and then drops down 
to ambient pressure at the outer edge of the contact zone . The 
mean pressure, also shown in Figure 17c, has been found by dividing 
the load supported by the trunk by the nominal contact area. The 
trunk load contribution was obtained by taking the total load and 
siabtracting the cushion load component found from the product of 
the measured cushion pressure and area. The actual trunk pressure 
is also plotted, and, as can be seen, it agrees closely with the 
mean pressure obtained from the force balance. As a final check, 
the load supported by the trvink and cushion as a function of the 
cushion-to-trunk pressure ratio is shown in Figure 19. The results 
coirparing the Hybrid Trunk Model with test data show that the theory 
predicts the actual loads and load sharing between the trunk and 
cushion quite accurately (within about 10 percent) , so that the 
contact pressure distribution model does not need further refine- 
ment at this time. 
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(b) Trunk Load Sharing 

Figure 19. Trunk Model Evaluation — Total Load and Load Sharing 
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Trunk Damping. - In dynamic operation, the trunk is deformed 
cyclically both in tension and flexure, and energy dissipation in 
the trunk material gives rise to a damping force which opposes the 
strain rate. Because the present trunk analysis does not solve 
for strain (and hence strain rate) , a damping model that links 
trunk material properties directly to trunk damping forces cannot 
be developed. An alternate approach, in which the damping charac- 
teristics are modeled by dimensional analysis (similarity) based 
on test data thus appears more appropriate. In keeping with the 
method of approach outlined earlier, the trunk is divided into 
segments (Figure 20) and a series of dashpots — one for each seg- 
ment — is included in the model such that the segment damping 
force is proportional to the vertical velocity of the 

trunk segment. 

Each dashpot B models the energy dissipation characteris- 
0 

tic of the trunk segment. Although all parts of the trunk dissi- 
pate energy, the major contributions will come from those parts 
that undergo high stress reversals, since the strain rate is high- 
est in these sections. Observations of a trunk in dynamic opera- 
tion suggest that the high stress reversal regions lie along the 
periphery of the trunk-ground contact zone, because it is here that 
the rate of change of trunk slope (and hence stress) is high and 
constantly changing with time as the contact area changes. As, 
a first order approximation, the damping model derived here 
assumes that all the energy dissipation in the trunk is concentra- 
ted along the trunk-ground contact periphery so that the damping 
coefficient of each dashpot B^ depends on the perimeter of the 
ground contact zone. This means that when a segment is not con- 
tacting the ground it has zero damping, and when it is contacting 
the ground it has a damping coefficient proportional to the contact 
perimeter. Thus, the damping coefficient of the dashpot associated 
with the i^ trunk segment is 
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Figure 20. The Trunk Damping Model 



( 30 ) 


B^(i) = DZ(i) 
e c c 

where D is a characteristic damping constant for the 

o 

trunk , 

and Z (i) is the perimeter of the ground contact zone 

for the i^ segment. 

The value of the damping constant will depend primarily 

on the material properties of the trunk, and correlation of exist- 
ing test data will provide an estimate of its value for the vari- 
ous types of trunks of interest. Although, at present, only 
limited data is available, the values of thus far obtained 

(see subsequent section on Parameter Evaluation) provide an ini- 
tial guide to the selection of this parameter. 

The damping forces and torques required to evaluate Eqns (22) , 
(23) and (25) are given by 


F 


ct 


E B (i) V. (i) 
. e t 


(31) 


T 


tx 


^ ®e \ 

1 


(32) 


T 


tz 


E B (i) V. (i) Z (i) 
. e t cx ' 


(33) 


where X and Z are the moment arms in pitch and roll, and, 

the minus sign is included because damping forces and 
torques oppose the velocity. 
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Model Synopsis 


The Flow System (Figure 7) 

(a) The fan is characterized by a static pressure rise 
element for forward and back flow in series with an 
inertance (duct) and a capacitance (volume) . 

(b) The trunk and cushion volume are found from the Hybrid 
Trunk Model, which characterizes the side trunk segment 
as an ideal two-dimensional membrane and the end 
segment as a "frozen” trunk. 

(c) The orifice areas between the trunk and cushion, trunk 
and atmosphere and cushion and atmosphere are found 
from the trunk shape as predicted by the Hybrid Trunk 
Model, along with the cushion orientation and ground 
profile . 

(d) The pressure within the cushion, trunk and plenum is 
considered to be uniform. 

(e) The pressure in the trunk/ground contact zone is found 
from the triangular profile given by the Hybrid Trunk 
Model . 

(f) The flow through the plenum, trunk and cushion is gov- 
erned by the unsteady state flow continuity equation 
in which the air is assumed to behave like a perfect 
gas and follow a polytropic expansion relationship. 

(g) The flow through all orifices is found from the incom- 
pressible flow square-law orifice equation. 

(h) The pressure relief valve is modeled as a second order 
mass-spring-dashpot system with nonlinear stops. 
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The Force System (Figure 8) 


(a) The mean contact pressure in the trunk/ground contact 
zone is equal to the trunk pressure. 

(b) The trunk contact area and location relative to the 
aircraft CG is found from the trunk shape predicted 
by the Hybrid Trunk Model. 

(c) The cushion area and location relative to the aircraft 
CG is found from the Hybrid Trunk Model. In width, the 
cushion extends between the lowest (ground tangent) 
points of the side trunk segments. In length, it ex- 
tends between the ground tangent points of the end trunk 
segments, or, if in ground contact, between the inner 
edges of the contact zone. 

(d) The total forces and moments acting on the aircraft 
occur due to the mean trunk contact pressure acting over 
the contact area, the cushion pressure acting over the 
cushion area, aerodynamic drag and trunk damping losses 
caused by aircraft heave motion, and trunk-ground fric- 
tion. 

(e) The forces and moments are found by dividing the cushion 
(and trunk) into segments, approximating the actual 
ground profile underneath the cushion by a similar set 
of segments parallel to the cushion, computing the 
cushion and contact pressure forces and moments for 
each segment, and then slimming them to determine the 
total force and moment about the aircraft CG. 

(f) The heave motion of the aircraft is found by applying 
Newton's law in the vertical direction to the aircraft 
CG. 
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(g) Angular accelerations in pitch and roll are obtained by 
applying the torque-angular momentum relation about the 
aircraft pitch and roll axes. 

(h) A coordinate transformation is carried out to express 
vehicle frame angular velocities in terms of Euler 
angles and their derivatives. 

(i) The angular velocities expressed in terms of Euler 
angles are integrated to give the angular position of 
the aircraft as a function of time. 

TESTING 


Test Objectives and Plan 


The object of the testing was to obtain data on the operation 
of a prototype cushion to help develop and verify the analytical 
model. To minimize hardware costs, it was decided at the outset 
to use an existing cushion available from NASA rather than build a 
new one. The test prototype selected was a unit built by Boeing 
under an earlier Navy contract. The main reasons for choosing this 
cushion for detailed study were that it was of a size suitable for 
laboratory testing, and had a very uniform molded plastic trunk 
free of seams, stitching or other anomalies. 

A test plan was drawn up with the following objectives in 

mind. 


1 . Parameter Evaluation 

• Trunk orifice discharge coefficient 

• Trunk damping constant 
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2 . Subsystem Model Development 

• Trunk shape in the side and end segments 

• Mean trunk-ground contact pressure and pressure 
distribution 

3 . Overall Model Verification 

• Static heave, pitch and roll stiffness 

• Peak values of pressures, displacement and accelera- 
tion in dynamic operation 

* 

The test plan is summarized in Table III. The tests in- 
cluded static tests, heave drop tests and pitch and roll moment 
release tests. The static tests provided key steady-state data, 
such as contact pressure distribution and trunk shape, needed to 
develop the trunk model. The dynamic tests provided time histories 
and peak values of pressures, displacement and acceleration which 
were used to verify the overall ACLS analysis. A key feature of 
the tests is that the modes are excited independently, so that 
coupling effects do not obscure the basic natural mode character- 
istics. A test rig was designed to carry out this plan by locking 
out various modes and thus allowing the cushion to be tested in 
pure heave, pitch and roll independently. A view of the test rig 
is shown in Figure 21. Details of the test hardware are given in 
the section that follows. 


The test parameters are actual test values. 
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Figure 21, The Test Rig 
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Test Hardware 


The Prototype Cushion 

The prototype ACLS selected for testing was a unit developed 
earlier by Boeing as part of a Navy contract*. A diagram of the 
ACLS including details of the trunk and fan are shown in Figure 22. 
The unit consists of an aliiminum frame which forms the plenum and 
provides the attachment points for the trunk. The fan is mounted 
directly on the plenum at one end. The trunk, which is made of 
fabric reinforced molded polyurethane, is attached at the inner and 
outer periphery by means of metal retaining strips bolted to the 
plenum, and was designed for easy removal and replacement. This 
turned out to be a very useful feature since the original trunk 
appeared to have been weakened through prolonged storage in the 
creased condition and ruptured repeatedly during the static pitch 
tests. The initial failures were repaired successfully using a 
commercially available two part polyurethane casting resin. Later, 
however, a major rupture developed in the inner attachment zone and 
the trunk had to be replaced. Fortunately, a back-up trunk was 
available. Static tests repeated with the second trunk gave results 
virtually identical to those obtained with the first trunk. 

Apparently due to better storage, the backup trunk survived all 
static and dynamic tests without rupture. 

The ACLS (with fan) had a mass of about 61 kg (135 Ibm.) . 
However, during testing a counterweight was added to balance the 
off-center weight of the fan and set the nominal pitch moment to 
zero. This counterweight along with the constraint mechanism of 
the test fixture increased the total mass to 89 kg (195 Ibm.). The 
rotational inertias of the system about the CG, calculated from 
measurements of its angular acceleration as a compound pendulum, were 


"k 

See Table II for reference. 
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(a) Overall View 


Figure 22. The Test Cushion 


(4.0) 





Approx. 1600 Holes 

Rib 

Detail I Section A-A 

Dimensions in cm (in) 


(b) The Trunk 


Figure 22 (continued) . The Test Cushion 
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2 2 2 
found to be 15.4 kg-m (11.38 slug-ft. ) in pitch and 4.53 kg-m 

2 

(3.35 slug-ft ) in roll. The fan, which requires a 200V, 400 Hz, 
3(j), 35 kVA power supply was powered by a motor- generator set. 

The Test Rig 

The test rig was designed to carry out static and dynamic 
testing of the ACLS both in the individual isolated modes (heave, 
pitch or roll) and in the coupled modes (heave-pitch or heave- 
roll) . An overall view of the test rig is shown in Figure 23. 

It consists of a stand made of Unistrut channels reinforced by 
cross members, and supporting an aluminum frame to which is at- 
tached a transparent plastic baseboard at eye-level. A 
non-rotating linear bearing is fastened to the frame at the 
center of the baseboard, and a steel shaft (heave pole) which 
slides through the bearing is attached at its upper end to the 
cushion via a pivot and yoke assembly. The system was designed 
so that the pivot axis passes through the CG of the cushion. The 
function of the yoke is to allow the natural modes of motion to 
be isolated as follows: 

(a) Pure Heave. When the lower part of the yoke is 

bolted directly to the top plate of the cushion (Fig- 
ure 24a) , the pivot is constrained from rotating, and 
the cushion can only move in heave as the heave pole 
slides through the linear bearing. Any desired pitch 
or roll angle for the cushion can be set by introduc- 
ing an appropriate shim or wedge between the yoke and 
top plate. In this mode, the cushion is raised or 
lowered by a hydraulic jack acting on the lower end 
of the heave pole. 

(b) Pure Roll. When the bolts fastening the lower yoke to 
the top plate of the cushion are removed, the cushion 
is free to rotate in roll about the pivot (Figure 24b) . 
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Accelerometer 



<T) 


(a) Side View (Section) 
Figure 23 . The Test Rig 



Pivot 



Figure 23 (concluded) . The Test Rig 





(b) Pure Roll 

Figure 24 The Yoke Assembly 


71 




In this mode, a collar clamped onto the heave pole 
just above the bearing sets the height of the cushion 
and prevents vertical motion. 

(c) Pure Pitch . When the yoke assembly is unfastened, 
rotated 90°, and then reattached, the pivot orienta- 
tion then allows the cushion to rotate freely in 
pitch. 

(d) Coupled Motion. Coupled heave-pitch or heave-roll 
motion is achieved by setting up the system in the 
pitch or roll mode and removing the collar clamped 
onto the heave pole to allow it to slide through the 
bearing. 

The Instriimentation 

The instrumentation used for the static and dynamic tests is 
summarized in Table IV and discussed below. 

Static Test Instrumentation. - In the static tests, pressures 
were measured with conventional U-tube manometers. In addition to 
measuring fan inlet and outlet pressure, plenum, trunk and cushion 
pressure, base pressure at five locations in the side and end 
trunk contact zone were also measured (see Figure 23b) . Displace- 
ments were measured with a scale: directly for linear motion and 

indirectly (by measuring the excursion of a radial arm) for angu- 
lar motion. The load on the cushion (including its self-weight) 
was found by allowing the lower end of the heave pole to rest on a 
spring scale. The trunk shape was determined by measuring the 
height of the trunk above the baseboard at discrete points. This 
was done by means of the eleven profile pins (knitting needles) 
spaced 3.2 cm (1.25") apart along the baseboard under the trunk as 
shown in Figure 23a. Details of the profile measurement scheme are 
shown in Figure 25. The pins can be moved up or down to touch the 
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TABLE IV 


SUMMARY OF INSTRUMENTATION 


Sensor No. 

Variable Measured 

Static Instrumentation 

Dyn£unic Instrumentation 





1 

Trunk Shape 

Profile Pins 

— 

2 

Trunk Contact Pressure 

Manometer 

— 

3 

Fan Inlet Pressure 

Manometer 

Pressure Transducer 

■ 4 

Fan Exit Pressure 

Manometer 

— 

5 

Plenuirii Pressure 

Manometer 

Pressure Transducer 

6 

Trunk Pressure 

Manometer 

Pressure Transducer 

7 

Cushion Pressure 

Manometer 

Pressure Transducer 

8 

Heave Displacement 

Scale (Ruler) 

Linear Potentiometer 

9 

Pitch Angle 

Scale (Ruler) 

Rotary Potentiometer 

10 

Roll Aiigle 

Scale (Ruler) 

Rotary Potentiometer 

11 

Cushion Load 

Spring Scale 

— 

12 

Heave Acceleration 

— 

Accelerometer 

13 

Pitch Acceleration 

— 

Accelerometer 

14 

Roll Acceleration 

— 

Accelerometer 


to 











I 


trunk, and are held in place by friction grommets. The positions 
of the lower end (head) of the pins, which are equivalent to the 
trunk height coordinates, are read off on a scale. The profile 
pin assembly is detachable, and can be positioned either under the 
side or end trunk segments. 

Dynamic Test Instriomentation . - A schematic diagram of the 
dynamic instrumentation system is shown in Figure 26. Signals 
from the transducers pass through signal processors and are re- 
corded on an eight channel strip chart recorder. Each signal 
processor consists of an amplification stage followed by a low-pass 
Butterworth filter to cut out noise (Figure 27) . Two types of 
processors were used depending on the strength of the transducer 
output signal. The pressure transducers and accelerometer, which 
have a low output, were connected to high gain processors (gain of 
100). The linear and rotary potentiometers, which have higher 
outputs, were used with low gain processors (gain from 1 to 10). 


The instrumentation specifications are given in Table V. 

Strain gage pressure transducers were used to measure system pres- 
sures and communicated with the plenum, trunk, cushion, etc. 
through short plastic tubes. Acceleration was measured with a 
moving mass linear accelerometer. To measure angular acceleration, 
this unit was mounted a known distance away from the center of 
rotation so that the output could be converted to angular accel- 
eration by dividing by the moment arm. The linear potentiometer, 
which measures the heave displacement of the cushion was attached 
vertically to the aluminium frame, with the moving arm (wiper) 
attached to the lower end of the heave pole (see Figure 23 ) • The 
rotary potentiometer, used to measure roll and pitch motion, was 
attached to a bracket on the yoke and connected to an extension of 
the pivot (see Figure 24b) . 
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Power Pickoff 


220V, 60Hz 


Figure 26 . Instrumentation System Schematic (Dynamic Tests) 









DC Gain 

R=270 Kfl 100 

R=1.8 KS2 1-10 

Max. Input Signal 

R=270 Kfi +0.1V 

R=1.8 +1V 

Supply Voltage +10VDC 

Input Impedence 27 

Roll-off Frequency 10 Hz 


Figure 27. The Signal Processor 
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SENSOR 


VARIABLE 

MEASURED 


MANUFACTURER 


SPECIFICATIONS 






Pressure Transducers 

Fan Inlet, 

Statham Instruments , 

PM 131 TC± 

(Four) 

Plenum, 
Trunk & 
Cushion 
Pressure 

Oxnard, California 

Range : + 

(-+ 

Accelerometer 

Heave , 

Statham Instruments , 

G312, A5a- 


Pitch and 
Roll Accel. 

Oxnard, California 

Range : ± 5 g 

Linear Potentiometer 

Heave 

Waters Manufacturing, 

LF-24 


Motion 

Wayland, Massachusetts 

Range : 0 . 

Rotary Potentiometer 

Pitch & 
Roll Angle 

Waters Manufacturing, 
Wayland, Massachusetts 



WP 7/8 M 


2.5-350 
17.23 kPa 

2.5 psi) 
5-350 


61 m (24") 









All signal leads were shielded. However, the presence of 
high-current, high-frequency power cables going to the fan made 
it apparent that a low-pass filter stage was needed before the 
signals could be recorded. Trunk flutter also induced noise in 
the signals. For these reasons, a Butterworth filter was includ- 
ed in the signal processor. The frequency response of the filter 
is shown in Figure 28. Its roll-off frequency (3 db bandwidth) 
was set at 10 Hz, since this is about the upper limit of the 
dominant ACLS frequencies of interest. 

Test Procedure 

The procedure followed to carry out the static and dynamic 
tests of Table III is given below. 

Static Tests 

(a) Adjust yoke to set appropriate mode (heave, pitch or 
roll) , and angle of attack (heave test only) . 

(b) Raise cushion with hydraulic jack under heave pole. 

(c) Turn on power to fan 

(d) Inspect trunk 

(e) Lower cushion slowly onto baseboard 

(f) Set CG height and lock heave collar (pitch and roll 
tests only) 

(g) Set cushion load (and moment) by adding (offset) weight. 

(h) Take trunk profile readings (side trunk) (Sensor #1) 

(i) Take base pressure readings under side trunk (Sensor #2) 

(j) Move profile pins and base pressure taps to end trunk 
location 
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(k) Take trunk profile readings (end trunk) (Sensor #1) 

(l) Take base pressure readings under end trunk (Sensor #2) 

(m) Read fan inlet and outlet pressure and plenum, trunk 
and cushion pressure (Sensor #s 3-7) 

(n) Measure vertical position or angle (Sensor #8, 9 or 10) 

(o) Remove cushion load 

(p) Unlock heave collar (pitch and roll tests only) 

(q) Raise cushion 

(r) Inspect trunk 

(s) Turn off power to fan. 

Dynamic Tests 

(a) Adjust yoke to set appropriate mode and angle of 
attack 

(b) Raise cushion to drop height with hydraulic jack 
(heave test) or set CG height and lock heave collar 
(pitch and roll test) 

(c) Turn on power to fan 

(d) Inspect trunk 

(e) Activate instrumentation and recording system 

(f ) Add unbalanced weight to cushion (pitch and roll tests 
only) 

(g) Release jack (heave test) or unbalanced weight (pitch 
and roll test) and allow oscillations to die down 

(h) Turn off instrumentation system 

(i) Unlock heave collar (pitch and roll tests only) 
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Raise cushion 


(j) 

(k) Inspect trunk 

(l) Turn off power to fan 

Test Results 
Parameter Evaluation 


Discharge Coefficient . - The discharge coefficient for the 
trunk holes can be easily found from the fan characteristic and 
the measured plenum and trunk pressures. When the cushion is 
high up off the ground (outside ground effect - OGE) , all the trunk 
holes discharge to the atmosphere. From the orifice equation, the 
discharge coefficient is given by 



where 


and 


Q 

A 




is the total flow 

is the total trunk hole area 

is the air density 

is the pressure drop across the trunk 
orifices 


The flow Q which is a function of the plenum pressure can 
be estimated from the static fan characteristic shown in Figure 10. 
Solving Eqn (34) for the high up (OGE) test condition gives a trunk 
hole discharge coefficient of 0.76. 

The discharge coefficient for cushion outflow through the 
clearance gap, however, cannot be found from the tests, because the 
gap area is unknown. However, since the trunk profile in the 
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cushion gap region is very similar to a two-dimensional converg- 
ing nozzle, the discharge coefficient for cushion outflow will 
be close to unity. Thus, when quantifying the analytical model, 
discharge coefficients of 0.76 and 1 have been used for the trunk 
and cushion gap orifices respectively. 

Trunk Damping Constant. - There are three damping mechanisms 
in ACLS: (a) backflow losses in the fan, (b) stagnation pressure 

losses due to cyclic flow through the orifices, and (c) energy 
losses due to cyclic deformation of the trunk*. When analyzing 
test data to evaluate the trunk damping constant (see 

Eqn. (30) ) , it is important to consider only those test modes 
where the first two components are negligible, so that the re- 
sults quantify the damping contribution of the trunk and not that 
of the total system. In heave vibration, all three damping com- 
ponents are present, because the fan can stall and the flow 
through the plenum and trunk orifices is changing. In pitch and 
roll, however, (as is shown later in Figure 36) the fan flow and 
plenum and trunk pressures remain constant, so that only the trunk 
damping component is present in these modes. Because the contact 
periphery in pitch is larger, it is more convenient to use the 
pitch test data to evaluate the trunk damping constant D^. This 
was carried out by solving Eqns (30) and (32) to match the cyclic 
attenuation in pitch angle response (Test No. 2.2, Table III), 

which gave a trunk damping constant D of 0.15 kPa*s (3.2 lb • 

2 ^ 
sec/ft ) . This value is for a fabric reinforced polyurethane 

trunk, which is non- stretching in tension, but very flexible in 

bending. Pitch data of another trunk (the Buffalo) were also 

analyzed to give a trunk damping constant of 8.4 kPa-s (176 lb • 

2 

sec/ft ) . This value is for a Nylon ply reinforced rubber trunk 


If present, relief valve flow losses will be a fourth damping 
mechanism. 
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made of two-way stretch material that is also flexible in bending. 
Although additional trunk data is needed before the form of the 
damping model can be verified, the two initial results show the 
right trend, because the more stretchable and larger scaled 
Buffalo trunk probably dissipates more internal energy and hence 
gives a higher value for the trunk damping constant. 


Trunk Shape and Base Pressure 

The test data confirmed that the trunk side segments behave 
differently from the end segments. As the cushion loading in- 
creases, the side segments bow outward while the end segments 
flatten (see Figure 12 ) . The base pressure distribution under 
the segments is also different. For the side segment, the cushion 
pressure drops down to ambient pressure very close to the ground 
tangent point of the trunk, as shown in Figure 29 . Thus, as men- 
tioned earlier, the ground tangent points form the side edges of 
the cushion. For the end segments, however, the base pressure 
under the trunk first increases and then drops down to the ambient 
level (Figure 30 ) . There is thus a ground contact region under 
the end segments which provides an additional force component. 

This type of trunk behavior has already been discussed, and, in 
fact, formed the basis for the Hybrid Trunk Model developed ear- 
lier. 


Static Characteristics 


The static characteristics of the cushion are shown in 
Figures 31 to 33 . Figure 31 shows the variation in cushion and 
trunk pressure as the load is increased. The cushion pressure. 


Unlike the trunk tested in this program, the Buffalo trunk is 
designed as an "elastic" trunk, which draws itself tight against 
the aircraft fuselage when not inflated. 
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Base Pressure Profile Coordinate, cm Base Pressure Profile Coordinate, cm 

kPa, gage kPa, gage 



(a) Load 867 N (195 Ibf) 



Distance from Cushion Center, cm (inch) 
(b) Load 1900N (427 Ibf) 


Figure 29. Trunk Profile and Base Pressure (Trunk Side) in Heave 
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Test #1.1 (Heave) 





which is zero at no load, increases almost linearly up to the 
maximiim test load. This increase in cushion pressure occurs because 
the cushion-to-atmosphere gap area reduces as the load is increased, 
so that the cushion pressure must rise to expel the inflow coming 
from the trunk. The trunk pressure also rises but begins to level 
off at the higher loads. This happens because the fan is now opera- 
ting near its maximTim pressure point (see Figure 10) and large flow 
changes can take place without much change in pressure. It must be 
noted that although the cushion-to-trunk pressure ratio increases 
with load, the portion of the load carried by the trunk remains fairly 
constant at about 30 percent (see Figure 19) . This happens because 
the trunk contact area increases rapidly with load, while the cushion 
area is substantially constant, so that the products of pressure and 
area increase at roughly the same time. 

In all tests, the plenum pressure was found to be equal to the 

trunk pressure because the plenum-to-trunk orifice is very large 
2 9 

(0.347 m compared to 0.0124 m for the trunk-to-atmosphere orifice 
area in OGE) . 

Figure 32 shows the hard surface clearance as a function of the 
load. In this case, the hard surface clearance is the distance 
between the trunk attachment points and the ground. The clearance 
decreases as the load increases due to compression of the trunk. The 
stiffness of the cushion, which is the inverse slope of the curve, 
starts initially at a low value and then increases as the trunk begins 
to stiffen. 

Figure 33 shows the pitch and roll angle as a function of the 
applied moment. The tests were carried out first by increasing the 
moment up to its maximum value and then reducing it. The data shows 
a little hysteresis, probably due to nonuniform retraction of the 
trunk caused by friction in the contact zone. As expected, the 
cushion is much stiffer in pitch than in roll, primarily because the 
restoring forces in pitch have a much longer moment arm, because the 
cushion is four times longer than it is wide. 
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Figure 33, Angular Deflection ^ 







Dynamic Characteristics 


The dynamic tests consisted essentially of heave drop tests 
and pitch and roll moment release tests. The chart recorder out- 
put for the 6" heave drop test (6" is the initial trunk-ground 
clearance) is shown in Figure 34. The first observation on this 
data is that the first four cycles following the drop are different 
from the remainder of the vibration* . This is because the initial 
impacts with the ground have sufficiently high kinetic energy to 
cause the fan to stall and delay the pressure rise. Each stall 
cycle, however, dissipates a part of the impact energy, so that 
after a few cycles, the remaining kinetic energy is below the fan 
stall threshold, so a non-stall vibration with faster pressure re- 
covery takes place. Another conclusion reached from the data is 
that heave damping is very low and after the major perturbations 
have attenuated, the system exhibits a low amplitude limit cycle 
instability . 

The stall behavior of the system can be explained with the 
help of Figure 35 . In Figure 35a , the operating point of the fan 
starts from A on the static characteristic and moves along locus 
A B C as the trunk hits the ground and the flow is reduced. When 
the trunk starts moving upwards, after the impact, the flow 
increases and the fan operating point moves along locus C D E. 

The cushion has now bounced back up, and the flow returns to its 
initial value, A. Since the pressure drop between the fan outlet 
and trunk is small, the trunk pressure is practically equal to the 
fan outlet pressure. The trunk pressure variation can then be 
plotted as shown in Figure 3 5b , which is very similar to the ob- 
served trunk pressure variation shown in the insert. Thus the 
first trunk pressure peak B corresponds to the maximum pressure 
the fan can generate when the flow is reducing, and the second 


The fifth cycle marks the transition between stall and non-stall 
vibration. 
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Fan Outlet Pressure 

Trunk Pressure (Trunk Pressure) 
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Fan Flow 

a) Fan Operation During Impact 



b) Trunk Pressure 

Figure 35. The Stall Behavior of the Cushion 


peak D corresponds to the maximum fan pressure when the flow is 
increasing. The fact that these two peaks have different magni- 
tudes is because, in one case, part of the static fan pressure 
rise is taken up in accelerating the flow, while in the other, 
deceleration of the flow adds an extra component to the static 
pressure rise. 

The trough C in the trunk pressure corresponds to the low 
flow point of the fan locus. This is verified by observing that 
the fan inlet pressure, which is a measure of the flow, has 
dropped sharply when the plenum pressure is at a minimum. 

This behavior seen in the drop tests confirms the earlier 
fan test results that formed the basis for the dynamic fan model. 

The stall behavior of the fan also causes the cushion pres- 
sure and acceleration during impact to have two peaks rather than 
one. As the cushion approaches the ground, the pressure and 
accelerating force begin to rise. Then, as the fan stalls, 

(locus B C, Figure 35a) the flow through the system reduces and 
the cushion pressure and accelerating force reduce. At point C 
on the fan curve, the pressure and acceleration reach a minimum, 
and then begin to rise as the fan recovers (locus C D) . Finally, 
when the cushion bounces back up, the gap area increases and the 
cushion vents to the atmosphere, causing the cushion pressure to 
drop to zero, so that the system moves under gravity acceleration 
only. The heave displacement, however, does not show a double 
peak during stalled ground impact because the force variation 
takes place too quickly for the system to respond. There is how- 
ever, a difference between the first four stall cycles and the 
remainder of the vibration; the difference being in the energy 
lost by damping. During each stall cycle, the system dissipates 
about 35 percent of its initial energy, as seen by comparing the 
successive heights of the rebound. This energy is dissipated 
through three mechanisms: hysteresis losses during fan stall, 

energy dissipation in the orifices, and flexural losses in the trunk. 
After about 80 percent of the drop energy is dissipated (in 
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about four cycles) fan stall does not occur, so one of the energy 
loss mechanisms is eliminated and system damping is reduced. 

In fact, when fan stall is eliminated, the net damping is so low, 
that the system is only marginally stable, and can be excited to 
execute low amplitude limit cycle oscillations, as seen in the 
latter part of the drop test record. 

Typical pitch and roll moment release tests are shown in 
Figure 36. The fan flow, as measured by the fan inlet pressure, 
remains constant, independent of angle, and so the plenum and 
trunk pressures also do not change. The cushion pressure, which 
starts out at zero because the cushion vents to the atmosphere 
when the pitch or roll angle exceeds a certain value, rises up 
to its equilibrium value as the moment is released and the cushion 
aligns itself parallel to the ground. The angular acceleration 
and displacement show a conventional damped vibratory response. 
System damping in pitch or roll is higher than in heave and a 
limit cycle instability does not occur. 


VERIFICATION OF AI-IALYSIS 

Computer simulations of the principal test cases of Table III 
were carried out to check the validity of the analysis and identi- 
fy the strengths and weaknesses of the analytical model. The 
results also give a rough idea of the accuracy that can be expec- 
ted from the analytical predictions. 

The Static Model 

Figures 37 and 38 compare theory and experiment for the 
static loading of the system in heave. The first figure shows 
the variation in hard surface clearance, and the second figure 
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Figure 36. Dynamic Test Record (Pitch and Roll) 
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Figure 37. Static Model Verification 
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shows the corresponding cushion and trunk pressures as the load 
is increased. The theory agrees reasonably well with experiment; 
the deflection and cushion pressure being within 10-15 percent 
from their test values, and trunk pressure predictions being even 
more accurate. Figure 39 shows the static pitch and roll angles 
as a function of the applied moment. The agreement here is not 
as good as in heave, with the analysis underestimating the angular 
deflection by 30-40 percent. The reason for this difference lies 
in the simplifications made in formulating the trunk model, in 
which the end segment is modeled as a frozen membrane due to the 
restraining effect of hoop tension. In fact, hoop tension drops 
off when the cushion is rotated in pitch or roll, so that the 
frozen approximation is less realistic in these modes than in 
heave because the reduced hoop restraint makes the trunk less 
stiff in angular deflection. 

The Dynamic Model 

Figures 40 and 41 compare the dynamic simulation results 
with the test records for the 15 cm (6") heave drop test. A 
10 Hz low pass filter was included in the simulation to account 
for the test filter used to cut out noise. The first figure 
shows the dynamic trunk and cushion pressure history following 
the drop. As mentioned earlier, the first four vibration cycles 
have sufficient energy to bring about fan stall, which gives a 
characteristic dual peaked shape to the trunk and cushion pres- 
sure waveform. Then, after sufficient energy has been dissipated 
so that the fan backpressure does not exceed the stall threshold, 
the remainder of the vibration takes place without stall. The 
analysis predicts this behavior remarkably well, including the 
dual peaked waveform and the transition from stalled to non- 
stalled operation. The period of the stall cycle is also accu- 
rately predicted. Figure 41 compares theory and experiment for 
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Figure 41. Dynamic Model Verification — Heave 
Acceleration and Displacement 
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the heave acceleration and displacement. Here, too, the analysis 
agrees well with the test. The acceleration waveform is similar 
to that of the cushion pressure because a major part of the ac- 
celeration is caused by the cushion pressure acting over the 
cushion area. Thus, the acceleration waveform also exhibits the 
characteristic dual peaked stall cycles followed by the unstalled 
vibration. The heave displacement shows the reduction in damping 
between the initial and final cycles of the simulation. Like the 
experiment, the initial cycles are more damped because of the 
energy loss associated with fan stall, while the remainder of the 
vibration has much less damping and finally settles down to a low 
amplitude limit cycle oscillation. 

Figure 42 shows the comparison between theory and experiment 
for the 41.4 mrad (2.35°) pitch moment release test. The analytical 
results confirm that fan dynamics in pitch are insignificant be- 
cause the trunk pressure (and hence fan back pressure) remains 
substantially constant. The simulation is able to predict the 
cushion pressure and displacement history very well, although 
peak acceleration estimates are higher than their measured values, 
apparently due to the simplifications inherent in the trunk dam- 
ping model . 

Although the complete time history of vibration provides very 
useful data on ACLS operation, the real value of the model as a 
design tool lies in its ability to determine the peak levels of 
the critical cushion parameters during touchdown and slideout. 

For instance, the maximum trunk pressure and dynamic trunk de- 
flection form direct inputs to the layout and structural design 
of the trunk. Similarly, the impact value of heave acceleration 
and pitch rotation provide a very good initial indication of the 
smoothness and stability of the landing. It is therefore appro- 
priate to evaluate the model based on its ability to predict the 
peak parameters of interest. Figures 43 to 47 show peak values 
and other key parameters as a function of the drop height (heave) 
or release angle (pitch and roll) . Figure 43 shows the maximum 
impact deflection of the cushion (below equilibriiim) , which is 
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Figure 43. Maximum Reduction in Hard Surface Clearance 
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Hard Surface Clearance Reduction, in 



equivalent to the maximum reduction in hard surface* clearance. 

When compared with the normal (static) clearance, this parame- 
ter will show how close the aircraft comes to hard surface 
contact during landing. As can be seen, agreement between 
theory and experiment is very good. Figure 44 shows the peak 
trunk and cushion pressure caused by impact. The simulation 
shows the right trends, although it tends to underestimate the 
cushion pressure by 10 - 25 percent. Figure 45 shows the peak 

acceleration and frequency associated with this acceleration 
(i.e., the stall frequency) . The frequency is predicted quite 
closely — within a margin of about 20 percent. The predicted accel- 
eration tends to be somewhat above the measured value, although the 
actual difference may be smaller than that shown because com- 
pliance in the members and joints of the aluminum cushion body 
(between the load points and the accelerometer mount) can reduce 
the measured acceleration to a level below that of a truly rigid 
structure. Figures 46 and 47 show the peak angular accelerations 
and natural frequencies in roll and pitch. The natural frequen- 
cies agree quite well; the theory predicting pitch frequency 
within a margin of 25 percent and roll frequency even more ac- 
curately. Roll frequency at the lower release angles is not 
shown because the test data was so damped that a complete cycle 
of oscillation was not available to estimate this parameter. Peak 
angular accelerations are not predicted as well as their corres- 
ponding frequencies, suggesting that the trunk damping model may 
need improvement. Unfortunately, as mentioned earlier, a more 
sophisticated damping model cannot be developed until a more 
advanced trunk deformation model is derived. Since the present 
trunk model gives good results for most of the other parameters, 
a major increase in simulation complexity may not be justified. 

A summary of the prediction error bound for the various parame- 
ters discussed above is given in Table VI. The static model values 
represent the differences between theory and experiment for the 
steady-state condition. The dynamic model values (except fre- 
quency) represent the corresponding peak levels caused by impact 
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Table VI 

Model Verification Summary 


Parameter 

Prediction Error Margin 

Static Model 

Dynamic Model 

Clearance Reduction 

10 - 15% 

<5% 

Trunk Pressure 

<5% 

10 - 15% 

Cushion Pressure 

10 - 15% 

10 - 25% 

Pitch Angle (Acceleration) 

30 - 50% 

(>50%) 

Roll Angle (Acceleration) 

20 - 50% 

(40 - 50%) 

Heave Acceleration 

N/A 

10 - 30% 

Stall Frequency (Heave) 

N/A 

10 - 20% 

Pitch Frequency 

N/A 

25% 

Roll Frequency 

N/A 

<5% 
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(heave drop or moment release) . With the exception of angular 
acceleration, which is difficult to predict because the trunk 
damping mechanism is very complex, the other outputs of the 
simulation compare well with test data; agreeing in most cases 
within a margin of 10-25 percent. The simulation could be further 
refined should a particular application warrant it to be, but 
even as it stands, it should serve as a valuable design and eval- 
uation tool for improving ACLS performance . 

MODEL USE 

In this section, some typical simulations of landing and 
taxi dynamics are presented, to show the kinds of results that 
can be obtained when the analytical model is used to simulate a 
full-scale ACLS aircraft. Rather than choose an arbitary con- 
figuration, the actual Buffalo aircraft was selected for simula- 
tion. The input data were collected from published reports (ref. 

5) and through discussions with the XC-8A (Buffalo) Project Office 
at Wright-Patterson Air Force Base . Simulations were carried out 
to determine 

(a) The static characteristics of the ACLS 

(b) The touchdown and slideout dynamics for a 41 m/s 

(80 knots) , 85 mrad (5°) nose up landing with a sink 
rate of about 1 m/s (3.5 ft/sec) , and 

(c) The behavior of the system during taxi at 3 m/s 

(6 knots) over an ungraded runway with a 23 cm (9 in) 
bump . 

Although the Buffalo ACLS is powered by two independent air 
sources, the simulations were carried out assuming that the sources 
shared the load equally. Also, the effects of trunk elasticity 
(see footnote on page 84) were not taken into account. The equi- 
librium characteristics of the Buffalo obtained from the static 
model are shown in Table VII. The slight offset between the 
aircraft CG and geometric center of the cushion gives the aircraft 
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Table VII 


PREDICTED EQUILIBRIUM FOR BUFFALO ACLS 
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a slight nose-up inclination at rest. The theoretica-1 fan power 
is the pumping power that would be required by an ideal fan. The 
actual power will thus be higher than this value due to losses in 
the fan. The fan stall margin is the maximum percent increase in 
backpressure that the fan can experience without stalling. Figure 
48 shows the static deflection of the Buffalo as the load and mo- 
ment are changed. The main conclusion emerging from the static 
analysis is that the aircraft has very lov>r stiffness in roll — 
a conclusion that was confirmed through the tests, and which led 
to a v;ingtip skid retrofit to increase roll stability. The rea- 
son for the shape of the roll deflection curve can be understood 
from the trunk model. ^*7hen a roll moment is applied, the first 
effect is to reduce the air gap under one of the side trunk seg- 
ments, while increasing the gap under the other. In this part of 
the operating characteristics the cushion does not generate any 
appreciable restoring moment, as shown by the initial vertical 
part of the curve. Then, as the gap under the lower side trunk 
segment becomes very small, the segment begins to bow outwards 
and subsequently touch the ground, thus increasing the restoring 
moment. At even larger angles, the cushion pressure begins to 
drop because the higher side trunk segment is well off the ground, 
thereby venting the cushion. In this operating region, the roll 
stiffness begins to reduce again as shown by the curve. 

The landing dynamics of the system are shown in Figure 49. 

From an initial nose-up attitude, the touchdown pitch moment 
causes the aircraft to pitch fonv^ard to a maximum nose-down pitch 
angle of 45 mrad (2-1/2°) . The pitch disturbance then dies out dur- 
ing slideout and braking. However, when the aircraft comes to a 
stop (about 9 seconds after touchdown) , the braking force (and 
moment) disappear and the aircraft undergoes a second pitch-heave 
vibration, coming to rest with the final nose-up attitude as pre- 
dicted by the static model. From Figure 49, the maximum increase 
in trunk pressure during landing is found to be ( 420-360) /360 = 

17 percent. Since the trunk pressure is practically equal to the 
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fan backpressure, it can be concluded that fan stall will not 
occur because the stall margin of 22 percent is not exceeded. 

The taxi dynamics of the system are evaluated by considering 
the case shown in Figure 50a, in which the aircraft crosses a 
ramp-like irregularity. This type of irregularity provides both 
gradual excitation and a sudden change . The irregularity length 
was chosen slightly smaller than the cushion length so that, for 
some duration, the entire obstacle would be contained entirely 
within the cushion. The heave and pitch motion response is shown 
in Figure 50b. Point A marks the spot where the leading edge of 
the cushion first touches the ramp. This causes the nose to 
pitch up and the aircraft CG to rise. Because the excitation is 
gradual, the change is relatively smooth. At point B, the front 
of the cushion goes over the edge of the ramp, so there is a 
sudden change in the pitch angle as the nose begins to drop down. 
The CG of the system also drops as the cushion pressure is re- 
duced by the sudden gap area increase at the front. At point C, 
the trailing edge of the cushion begins to move up the ramp, 
thus raising the system CG but continuing to increase the nose- 
down pitch angle. Finally, at point D, the back of the cushion 
goes over the edge of the ramp thus causing a sudden drop in CG 
elevation and restoring the aircraft to its normal pitch orien- 
tation. When crossing the obstacle, the aircraft CG rises by a 
maximum of about one-half the obstacle height, and the peak (nose 
down) pitch angle reaches a value of about 45 mrad (2-1/2°) . 

The main conclusions that emerge from the Buffalo simula- 
tions are as follows. 

• The aircraft will have poor roll stability. 

• In equilibrium, the aircraft will maintain a slight 
nose-up attitude. 
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• The theoretical (minimum) fan power (both fans com- 
bined) will be about 800 kW (1000 hp) . 

• The maximum (nose-down) pitch angle during landing 
will be about 45 mrad (2-1/2°) • 

• Under normal landing conditions, fan stall will not 
occur. 

• Peak landing accelerations will be on the order of 

1 g 

• The trunk will have to withstand peak pressures of 
about 20 kPa (3 psi) . 

• Although taxi behavior depends on the field roughness, 
the satisfactory negotiation of a 23 cm (9 in.) bump 
suggests that the aircraft will be able to operate 
from rough fields. 

CONCLUSION 

A heave-pitch-roll analysis and computer simulation has been 
developed to evaluate the landing and taxi dynamics of an ACLS 
aircraft. Initial results have shown that fan dynamics and stall 
have a major effect on system performance. Comparison with zero 
forward speed test data shows that the model can predict most of 
the key system parameters within a margin of 10-25 percent. 

The simulation can now be used as an analytical tool, to 
evaluate existing configurations and develop new designs that 
overcome current shortcomings such as low angular stiffness and 
poor stability. Concurrently with design development, additional 
test data should be obtained in two important areas, to complete 
the verification of the analysis. These areas are: 
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The stall and backflow characteristics of ACLS fans 


• The effects of forward speed on ACLS dynamics 

When this work is complete, NASA will have available the 
technology base needed to implement practical and efficient ACLS 
for any aircraft application of interest. 
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APPENDIX 

THE EQUATIONS OF THE HEAVE-PITCH-ROLL MODEL 


Aj^ to 


af 

A 

c 

^cn 

^dd 

^df 

A 


pa 


pc 

^Pt 

A. 


ta 

^tc 

A 

u 


^11 


Symbol List 
A 

Horizontal distance between inner and outer 
trunk attachment points 

Area 

Area of fan inlet orifice 
Cushion area 

Area of trunk in ground contact 
Area of duct connecting fan to plenum 
Crossectional area of fan flow passages 
Orifice area, cushion to atmosphere (gap area) 
Area of trunk hole 

Orifice area, plenum to atmosphere 
Orifice area, plenum to cushion 
Projected heave area of aircraft 
Orifice area, plenum to trunk 
Trunk crossectional area 
Orifice area, trunk to atmosphere 
Orifice area, trunk to cushion 
Upstream orifice area 
Relief valve orifice area 
Trunk crossectional area components 
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B 

b - Vertical distance between inner and outer 
trunk attachment points 


^11 

^12 

^13 

^21 

^22 
etc . 


Euler angle transformation factors 


B - Damping coefficient of trunk segment 


C 


'af 

C_ 


'du 

'fx 

'fz 

C 

g 

1 

'pa 


pc 


'Pt 


'ta 


Fan inlet orifice discharge coefficient 

X axis distance between aircraft CG and cushion 
center 

Drag coefficient for heave motion of aircraft 

Discharge coefficient of upstream orifice 

X axis distance between aerodynamic drag center 
and aircraft CG 

Z-axis distance between aerodynamic drag center 
and aircraft CG 

Discharge coefficient for cushion to atmosphere 
flow 

Discharge coefficient for plenum to atmosphere 
flow 

Discharge coefficient for plenum to cushion flow 

Discharge coefficient for plenum to trunk flow 

Discharge coefficient for trunk to atmosphere 
flow 
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- Discharge coefficient for trunk to cushion flow 

- Relief valve discharge coefficient 

D 

d - Distance between inner trunk attachment points 
dx - Width of straight trunk segment 
D^ - Trunk damping constant 

F 

F - Force 

F - Equilibrium vertical cushion force 
cn ^ 

F - Force on aircraft due to cushion pressure 
cp ^ 

F^^ - Trunk damping force along vehicle y-axis 

F^y - Total force along vehicle y-axis 

F^^ - Aerodynamic drag force along vehicle y-axis 

F^ - Z-axis distance between aircraft CG and center 
of cushion 

F^ - Force on aircraft due to trunk contact pressure 
tp 

F^ - Force acting along inertial X-axis 

Fy - ACLS force in the inertial vertical direction 

F^ - Force acting along inertial Z-axis 

G 

g - Gravity acceleration 

G - y-axis distance between aircraft CG and cushion 
^ center 
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H 

H - Angular momentum vector about the CG 

H - Cushion width 
w 

H . - Cushion width at no load 

W1 

H - Angular momentiim component along vehicle x-axis 

X 

H - Vertical distance between hard surface and 
^ lowest point of trunk 

Hy - Angular momentum component along vehicle y-axis 

H . - Value of H at no load 
yi y 

- Angular momentum component along vehicle z-axis 


I 


[l] - Inertia matrix 

- Inertance of air in fan ducts 

I - Segment location identification number 

- Segment type number 

- Moment of inertia of aircraft about the vehicle 
x-axis 


I 

y 


Moment of inertia of aircraft about the vehicle 
y-axis 


I - Moment of xnertia of aircraft about the vehicle 
z 

z-axis 


I \ 

xy I 

^yz I Cross products of inertia of the aircraft with 
I respect to the x, y and z axes of the vehicle 

I 1 

zx / 
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K 


k - Polytropic exponent for air expansion 

- Relief valve spring stiffness 
K - Trunk height design constant 

L 

- Perimeter of trunk crossection 

H - Perimeter of the ground contact zone for each 
trunk segment 

5,^^ - Length of duct between fan and plenum 

- Length of fan flow passages 

& - Peripheral distance from inner trunk attachment 

^ point to first row of trunk holes 

- Peripheral length of trunk sector (cushion side) 

1 2 - Peripheral length of trunk sector (atmosphere side) 

- Distance of hard surface from ground datum 
Lg - Length of straight segment of cushion 

M 

m^ - Mass of pressure relief valve 

M - half the number of segments in each straight 
section of the trunk 

- Aircraft mass 
a 

N 

N - Half the number of segments in each curved 
section of the trunk 

- Number of trunk orifices per row 

- Number of rows of trunk orifices 
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N- 


V 




Q 



Q 

Q 

Q 

Q 

Q 

Q 

Q 


fp 

fx 

pa 

pc 

pt 

ta 

tc 

Qv 


Number of relief valves 

Number of trunk orifice rows communicating with 
cushion 

Number of trunk orifice rows communicating with 
atmosphere 


P 

Atmospheric pressure 

Fan inlet pressure 

Average contact pressure 

Cushion pressure 

Static pressure rise across fan 

Plenum pressure 

Relief valve preload pressure 

Trunk pressure 


Q 

Flow 

Flow from cushion to atmosphere 
Fan flow 

Flow from fan to plenum 

Dynamic fan flow 

Flow from plenum to atmosphere 

Flow from plenum to cushion 

Flow from plenum to trunk 

Flow from trunk to atmoshpere 

Flow from trunk to cushion 

Flow through pressure relief valve 



I 


R 

r - Radius 

R - Radius of equivalent circular crossection trunk 
R^ - Radius of trunk sector (cushion side) 

R2 - Radius of trunk sector (atmosphere side) 


S 

S - Peripheral length of cushion 
Sj^ - Spacing between the rows of trunk holes 
- Relief valve peripheral length 


t 

T 

T 


T 

T 

T 

T 


cpx 

cpz 

dfx 

df z 

T, 


' f z 


T. 


nx 


T 


nz 


T4. 

tpx 


T 


tpz 


T 


Time 

Hoop tension in trunk 

Torque vector about the aircraft CG 

Cushion pressure torque about vehicle x-axis 

Cushion pressure torque about vehicle z-axis 

Aerodynamic torque about vehicle x-axis 

Torque due to aerodynamic drag 

Torque due to trunk contact friction 

Cushion and contact pressure torque about vehicle 
x-axis 

Cushion and contact pressure torque about vehicle 
z-axis 

Trunk contact pressure torque about vehicle 
x-axis 

Trunk contact pressure torque about vehicle 
z-axis 
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Trunk damping torque about vehicle x-axis 
Trunk damping torque about vehicle z-axis 
Torque component about vehicle x-axis 
Torque component about vehicle y-axis 
Torque component about vehicle z-axis 

U 

Unit vector along intermediate axis during 
Euler angle rotation 

Unit vector along inertial Y-axis 

Unit vectors along vehicle x, y, and z axes 


V 

Heave velocity component of aircraft 
Cushion volume 
Plenum volume 
Trunk volume 

Vertical velocity of trunk segment 

Velocity of trunk (contact) segment center 
relative to aircraft CG 

X 

Vehicle roll axis 

Relief valve stop clearance 

Relief valve motion 



X - Inertial roll axis 

X - Aircraft CG motion in the inertial forward 
direction 

X , - Longitudinal distance between center of pressure 

° and cushion center 

X - Longitudinal distance between segment center and 
cushion center 


X - Horizontal distance from centroid of trunk 

crossectional area to inner trunk attachment 
point 


X 


g 


Ground profile location coordinates 


X^ - X-axis distance between the center of the trunk 
contact area segment and the center of the 
cushion 


X^j^ - Lateral distance between center of trunk contact 
pressure and cushion center 


Xl to 


X 


cr 


X, ,x ,x 

b e' er 


X coordinates of centroids of trunk cross- 
sectional area components 


Y 

y - Vehicle yaw axis 
Y - Inertial yaw axis 

Y - Aircraft CG displacement in the inertial vertical 

direction 

Yg - Ground profile height coordinate 

Y , - Hard surface clearance 

gh 
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z 


z - Vehicle pitch axis 

Z - Inertial pitch axis 

Z - Distance along inertial Z axis 
eg ^ 

Z^j^ - Lateral distance between center of cushion 
pressure and cushion center 

Z - Lateral distance between segment center and 
cushion center 

Z - Ground profile location coordinate 

g ^ 

Z^ - z-axis distance between the center of the trunk 
contact segment and the center of the cushion 

Z^j^ - Lateral distance between center of trunk contact 
pressure and cushion center 

Z^ - Relief valve damping constant 

Greek 

a - Fan curve polynomial coefficient 

3 - Angle subtended by curved trunk segment 

6 - Angle between cushion axis and curved trunk 
segment 

6 - End trunk excursion 

e 

(Sg - Side trunk excursion 
p - Air density 

0 - Roll angle about vehicle x-axis 
6^ - Eulerian roll angle 
y - ACLS braking coefficient 
(p - Pitch angle about vehicle z-axis 
(j)g - Eulerian pitch angle 
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\p - Yaw angle about vehicle y-axis 


ip^ - Eulerian yaw angle 
w - Angular velocity 



Angles subtended by trunk sectors 


Other 

( ) - Velocity 
( ) - Acceleration 

(i) - Refers to i^ segment of trunk 
("^) - Vector 

( - Refers to initial or equilibrium value 

( ) . - Refers to value of variable when trunk is not 
1 

in ground contact 

( ) - Refers to change in value of variable caused 

by trunk-ground contact 

- Differentiation with respect to time 
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THE STATIC MODEL 


The Trunk Model 

To facilitate the analysis, the trunk is divided into eight 
sections as shown in Figure 51. Each section is divided into 
segments: M segments per straight section and N segments per 

curved section. Thus the total number of trunk segments is 
4 (M+N) . The location of each trunk segment (i.e., right or left, 
front or rear) can thus be identified by the section number I . 

5 

The type of each segment (i.e., straight or curved) is denoted 
by the binary variable I^. A value of = 1 represents a curved 
segment, while = 0 represents a straight segment. Numbering 
each segement clockwise starting from the center of the rear end 


■unk gives 

the 

following 

values of I and I . . 

s t 

Section No. 

Type 

Segment 

NO. (i) 


(V 

0 

< 

i 

< 

N 

1 

1 

N 

< 

i 

< 

N + M 

2 

0 

N + M 

< 

i 

< 

N + 2M 

3 

0 

N + 2M 

< 

i 

< 

2N + 2M 

4 

1 

2N + 2M 

< 

i 

< 

3N + 2M 

5 

1 

3N + 2M 

< 

i 

< 

3N + 3M 

6 

0 

3N + 3M 

< 

i 

< 

3N + 4M 

7 

0 

3N + 4M 

< 

i 

< 

4N + 4M 

8 

1 


The above grouping is subsequently used to keep track of 
the location and type of the segments when determining the 
moments generated by the cushion. 
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Front 



Figure 51. Trunk Sections (Bottom View) 
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Crossectional Shape 


Without Ground Contact - (a) The Side Trunk - The shape of 
the side trunk, made up of two circular arcs (Figure 52a) , is 
defined by the four independent parameters a, b, t, and Hy. The 
first three parameters are constants for any given design, while 
the fourth is found from the pressure ratio through a membrane 
analysis (ref. 11) of the side trunk. The trunk height can 
be expressed as 

H /H . = f (P /P. ) (A-1) 

y^ yi c^ t 

where 

H = trunk height 

y ^ 


H . = trunk height at P = 0 , and 

yi ^ c ' 

P^/P^ = cushion-to-trunk pressure ratio 

The above functional relationship is found by simultaneous 
solution of the following ten equations; the first nine for 
geometric compatibility and the tenth for force equilibrium. 


(A-2) 

^2 ‘^2 ~ ^2 (A-3) 

= Si (A-4) 

Cos = (R 2 ~Hy)/R 2 (A-5) 
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B^Cos (<j)^-Tr/2) + R 2 Sin ((>2 = a 

R^Sin((J)^-iT/2) - (Hy-Rj^) = b 



(A-6) 

(A-7) 


(A-8) 


(A-9) 




(A-10) 


(b) The End Trunk - The shape of the end trunk is also made 
up of two circular arcs and defined by the four independent 
parameters a, b, £, and Hy. However, because the end trunk 
shape is independent of the pressure ratio (frozen) , the trunk 
height Hy is constant. 

H = K (A-11) 

y 

where K is the design value of the end trunk height. For the 
simulations described in this report, H is taken to be equal 
to Hy^, and the trunk parameters are determined from equations 
(A-2) to (A-9) . 
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with Ground Contact - When the value of Hy, as calculated 
from Equation (A-1) or (A-11) , exceeds the available ground 
clearance, the trunk shape changes due to ground contact. In 
keeping with the assumption of the hybrid trunk model, ground 
contact only affects the shape in the contact zone, and causes 
the trunk to conform with the ground contour as shown in 
Figure 52b. In this situation, by definition of H 


H 




(A-12) 


However, the parameters describing the undeformed part of the 
trunk are determined from the calculated value of H . 


Segment Center Distance 

The segment center distance X (i) and Z (i) (Figure 4) are 

OX ox 

found as follows. 

Section 1 (1=1) 

X^j^(i) = -Lg/2 - (d/2 + R 2 Sincf> 2 ) Cos6 (i) (A-13) 

Z^^(i) = (d/2 + R 2 Sin(|) 2 ) Sin 6(i) (A-14) 

where 


6(i) =. (i-0.5)3, and 

3 = tt/2N 

Section 2 (I = 2) 
s 

(A-15) 
(A-16) 


X^^(i) = -L^/2 + (i-0.5-N) dx 

s 

Zcx(i) = ^2 Sin<|)2 + d/2 
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where 


dx = Lg/2M 
Section 3 (1^ = 3) 


X (i) 
cx 

(i-N-M-0.5) dx 


(A-17) 

Z (i) 

cx 

R 2 Sin<j)2 + d/2 


(A-18) 

Section 4 (I^ 

= 4) 




Lg/2 + (d/2 + R 2 

Sin(|> 2 ) SinS (i) 

(A-19) 

= 

(d/2 + R 2 Sin(Ji2) 

Cos6 (i) 

(A-20) 

6 (i) = (i-N- 

2M-0. 5) e 



Section 5 (I 

s 

= 5) 



X (i) = 

cx 

Lg/2 + (d/2 + R 2 

Sin(fi 2 ) Cos6 (i) 

(A-21) 

Z (i) = 

cx 

-(d/2 + R 2 Sincj)2) 

Sin6 (i) 

(A-22) 


where 

6(i) = (i-2N-2M-0.5) 6 
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6 ) 


Section 6 (I = 
s 


X (i) 
cx 

Lg/2 - (i-3N-2M-0. 5) dx 

(A-23) 

Z (i) = 

cx 

-(R2Sin<|)2 + d/2) 

(A-24) 

Section 7 (1^ = 

7) 



-(i-3N-3M-0. 5) dx 

(A-25) 

= 

-(R 2 Sin())2 + d/2) 

(A-26) 

Section 8 (I = 
s 

8) 


X (i) 
cx 

- L /2 - (d/2 + R Sincj>_) Sin5(i) 

(A-27) 

Z (i) = 

cx 

-(d/2 + R 2 Sini|) 2 ) Cos6 (i) 

(A-28) 


where 

6(i) = (i-3N-4M-0. 5) e 

Hard Surface Clearance 

The hard surface clearance for each segment, Ygj^(i) , is 

found from the position of the ACLS (X , Y , e , * ) and from 

eg' eg' e' ^e 

the ground profile (Figure 53) . This is carried out in three 
steps. 

(a) The parameters X (i) , Z (i) , ^-^Ci) and the orientation 
of vector DA for each segment are calculated 
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(b) The ground profile coordinates and DA are used -to 
calculate the length CA 

(c) The hard surface clearance for each segment is then 
given by £^(i)-CA 

The results of the analysis give 

^ - = c ] >^11 

- [»4<i> + Gg] b2i + - pj Si 

"g<i> = ‘>13 

+ <=g] 1=23 ^ [2ox<^> -*■£] >=33 

* 4 '^* “ '^ og ^*^22 ■'■ " “^ o ] ^ 12^^22 

* C^cx'i) - ^f] ‘>32/t'22 - % 


where b^ the transformation factors that convert from the 
vehicle frame to the inertial (Euler Angle) frame are given by 


11 

COSllJ 

e 

CosiJ) + Sine Sin* Sin* 
e e ^e e 

12 

Sin<})^ 

Cose 

e 

13 

-Sin4) 

( 

Cos* + Sine Sin* Cos* 
s e ■ e e e 

21 

Sini^^ 

Cos4> Sine - Sin* Cos* 
e e ^e e 
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t>22 “ ^ Cos(()^ Cos6^ 


= Sincj) Sinil> + CosTp Cos(j> Sind 

O 6 6 6 6 


The length CA is given by 


CA = Y (i) Cosd) Cos6_ 
g ' ^e e 


CA-32) 


where the ground coordinate ^g(i) is found from the profile 
relationship 


(i) = f [x (i) , Z (i)] 


(A-33) 


The hard surface clearance is then given by 


Ygj^(i) = £^(i) - CA 


(A-34) 


Areas and Volumes 


The orifice areas and cushion and trunk volumes, for a 
particular trunk orientation, are calculated independently for 
each segment and then combined to give the total system value. 

For convenience, some of the areas and volumes are divided into 
two components - the i component and the r component (denoted by 
the subscripts i and r) . The i values are calculated assuming 
that the trunk segment under consideration is out of ground 
contact. The r values represent the changes in the segment areas 
and volumes due to trunk-ground contact. The actual segment 
areas and volumes are fo\ind by subtracting the respective r 
values from the i values. The total areas and volumes are deter- 
mined by combining the areas and volumes for each segment. For 
example 
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I 


V 


t 


4 (M+N) 

i=l ' ' 


CA-3 5) 


where 


V, 


total trunk volume 


v^i(i) 


V^r<i) 


(■ 




= i value of trunk volume for segment 

t h 

= r value of trunk volume for i ■ segment 

th \ 

= 0 if segment is not in ground contact.! 


Similar relations hold good for cushion area, ; cushion 
volume, V^; and gap area, A^. Other parameters, such as , 
A^g^, A^^^, A^^^, and A^^ are calculated directly in a single 
step. 


Without Ground Contact - From Figure 54, the trunk cross- 
sectional area A^^(i) is given by 


where 


Ati(i) = A^ - A^ + A3 - A^ + A5 


(A-36) 


^1 = 

'^2 ^ 2 ^^ 

^2 

(R 2 -Hy) 

A 3 = 

4.^rJ/2 

^4 = 

Xb/2 


143 
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Sector 


Sector 


1 

= ABD 

6 

= BDJ 

2 

= ABC 

7 

= BJL 

3 

= EOF 

8 

= DFK 

4 

5 

= EGH 
= CFG 

9 

= FLK 


Figure 54 (concluded) , Trunk Shape Model 







and 


= (a - R 2 , Sin<t >2 - X) (Hy-R^) /2 


X 


b (a 
b 


- R2 Sinc|)2) 



The trunk volume for the 


. th 
1 


segment is 


given by 


Vti(i) 


dxA^.(i) ; = 0 

(2 + ' ^t ~ ^ 


(A-37) 


where X is the horizontal distance of the centroid of the area 
e 

A^^ from the inner trunk attachment point. X^ is calculated as 
follows 


X 


AjXi - A2X2 + A3X3 - A^X^ + A3X3 


(A-38) 


where X^, X 2 , etc., are the X coordinates of the centroids of 


the areas A^ , A 2 f etc., respectively. 


Xn = R 


2 Sin (}.2 - 4Sin (<}> 2 / 2 ) 


X 2 = 0.6667 R 2 Sin(f )2 

X^ = R 2 Sin <})2 + 4Sin^ R^/3(j)^ 


X 4 = a - 0.333X 


X^ = R 2 Sin (|)2 + 0.333 (a - R 2 Sin 4>2 - X) 
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The cushion area for the segment is given by 




|(d/2 + R2Sinc|)2) dx ; = 0 

l(d/2 + R2Sin<})2)^ 6/2 ; 1^=1 


CA-39) 


To calculate the trunk-to-cushion flow area it is necessary 
to start with the number of trunk holes communicating with the 
cushion. The number of rows of holes communicating with the 
cushion is given by the integer value of ( 

number of communicating holes is /S, ) + 1 N, . The trunk- 

th ^ ^ P -in 

to-cushion flow area for the i segment is thus given by 


^tci(i) 


= integer 


-|-£ + 1 


^h^h 


dx/S (A-40) 


and 


for 1^=0 




- integer 


^2~^p 


+ 1 


^h^h 


6 (d/2 + R2Sin<})2)./S 


(A-40) 


for = 1 


where S, the cushion periphery, is given by 


S = 2Lg + 2Tr(.d/2 + R 2 Sin<J> 2 ) 
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The trunk-to-atmosphere flow area for the segment is 

given by 


^tai^i) = . dx/S - (i) 


(A-41) 


for = 0 


and 


^tai^^> = N^Vh 


P(d/2 + R2Sinct)2) 


-A^^i(i) (A-41) 


for = 1 


The cushion-to-atmosphere flow area (gap area) for the 

. th . 

X segment xs 


pgh<^> - «y] 


'• It = “ 


A . (i) = 

gi 


(A-42) 


[Ygj^(i) - hJ (d/2 + R2Sin(()2)B; ^ 


*th 

Finally, the cushion volume for the i^ segment is given 


by 


V . (i) 
cx 


Ygj^(i) (d/2 + R2Sin(()2)- + A 2 


dx (A-43) 


for = 0 


and 


= V‘i’ t RjSin+j)^ 


- B(d/2 + X^2^ ' ^t " 


(A-43) 
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where 


, _ “ ^2^2 

^12 - A2 


(A-44) 


With Ground Contact - Ground contact occurs when Y ,{i)< H 
for the trunk sides, or Ygj^(i) < K for the trunk ends. With 
ground contact, the trunk cross-sectional area changes as 
follows : 


A , ( i ) = 

tr 


Ag - A? + Ag Ag 


; ± tt/2 (A-4 5) 


Ag - A^ - A^g + A^^ ; 4>4 > tt/2 (A-46) 


where Ag , Ay, etc., are the areas of the sectors shown in 
Figure 54. 


* 3/2 


A_ 


(^2- "y ^h<Tl) 


R 2 ®i "*3 


RJ + 4/2 


(e, - H + Y ^(i)l 

Aj = 1 3ll_^^Sin+^ 

'10 = [Ri - «y + ■'l 


"11 


R^tt/4 
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and 


= Cos 


(j), = Cos 


-1 


-1 


‘2 - K - 



The r value of the trunk volume is calculated as 


follows 


r 


V . ( i) 

tr 


(Ag - + Ag - Ag) dx ; 

4)4 £ tt /2 
It = 0 


(Ag - A^ + 2Ag - 2Ag + A^q - A^^) dx ; 

<1>4 > tt/2 

It = 0 


i (Ag - A^ + Ag - Ag)B(X^^ + d/2) ; (A-47) 

<t>4 1 '^/2 
"t = 1 


(Ag - A^) 3 (2X^-^Xj^^ + d/2) 


+ 3d/2(2Ag - 2Ag + Aj^q - 




+ 2SX^^(A3 - Ag) - eX^^(Aj^j^ - Aj^g); 

■^4 ^ 


"t = 1 
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are the X-coordinates of the centroids of 


where and Xj^^ 

area 


er 


^ 6^6 ~ •^ 7^7 ^ ^ 8^8 ~ '^ 9^9 
^6 “ ^8 ” "^9 


^ 6^6 “ -^ 7^7 ^ 11^11 ^ 10^10 



br 

^6 

- ^ ^11 

0 
I — 1 

< 

1 


X^ , etc . , are X 

coordinates 

of the centroids 

etc . 

, respectively. 



^6 

= R 2 Sin<t )2 

- 4Sin^ (<t>^/2) 


^7 

= R 2 Sin <})2 

- 0 

. 333R2Sinc|)2 


^8 

= R 2 Sin (|)2 

+ 4Sin^ ( 4 > 4 / 2 ) 

Ri/3^4 

X 9 

= R 2 Sin <|)2 

+ 0 

. 333Rj^Sin<j)^ 


^10 

= R 2 Sin<f >2 

+ R 

1/2 


^11 

= R 2 Sincf )2 

+ 4Sin^(ir/4) R^ 

/1.5tt 

The 

r value of the 

cushion area 

A^^(i) is given 



/ 

dxR 2 Sin <()2 

; = 0 


A (i) = ^ 

cr ’ 

1 . 

1 ^ 

2 + R 2 Sinc|) 2 ^ R 2 Sin(J )2 



^ ■ 

2 2 
R 2 Sin ((>2 

; It = 1 


(A-48) 


(A-49) 
areas Ag , 


(A-50) 
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The mamber of trunk orifice rows communicating with the 
cushion is given by 


N 


1 


integer 



‘*’3^2 


+ 1 


(A-51) 


The area of the trunk orifices communicating with the cushion 
(for the i^ segment) is given by 


Atoi'i' 


NfNj^Ah dx/S 

NiNhAh3(d/2 + R2Sin(|)2)/S 


0 

(A-52) 

1 


Similarly, the number of orifice rows communicating with 
the atmosphere N 2 is given by 


N, 


- I + I + (N -1) S, - 4>.Rt 
1 p r ' h ^4 1 


+ 1 (A-53) 


The value of the trunk-to-atmosphere area for the i^^ segment is 


\ai 


^2^h^h 


; 1^ = 0 


(A-54) 


N2NhAh3(d/2 + R2Sinc[>2)/S ; ^ 


The orifice area communicating with the cushion for the i 
segment is given by 


. th 


^tcr^^^ = integer ( P + l) Nj^A^^ dx/S 


- -^tci^^^ ^t = ° 


(A-55) 
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(i) = N N,A, dx/S - A.^(i) - A. , (i) 

tar r h h tc ta 

- for It = 0 (A-56) 

\ar^^^ = N^Nj^Aj^a (d/2 + R2Sin<))2)/S - A^^(i) 

The r value of the clearance gap area is given by 

The above equation is easily derived when it is recognized that 
ground contact blocks off the gap area. 

The r value of the cushion volume is given by 

- dx (Ag-A^) ; = 0 

V^^(i) = (A-58) 

- a (d/2 + X^j.) (Ag-A^) ; = 1 
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where 


X, 


cr 


^ 6^6 ~ ^ 7^7 
Ag - A^ 


The trunk-ground contact area on the cushion side of the 
contact zone is given by 


A = R-Siniji^ dx for 1=0 

cm z ’ j ^ 


CA-59) 


A • 
cm 


1 

2 


(d/2 + R2Sin<))2) 


- (d/2 + R2Sin<()2 - R2Sin<|)2) 


(A-59) 


for = 1 

The trunk-ground contact area on the atmosphere side of 
the contact zone is 


A (i) = R,Siniti. dx 

cnr 1 ^4 


(A-60) 


for =0, <)>4 :i •^/2 




(d/2 + R2Sin(j>2 + R;j^Sin<J>^) 

- (d/2 + R2Sin(j)2)^J 


(A-60) 


for = If <(>4 5 ‘^2 
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A (i) = R, dx 

cnr ' 1 


(A-60) 


for = 0, > ir/2 


A ^(i) = 

cnr 


£ 

2 


(d/2 + R2Sin<|)2 + R;j^) 
2 


- (d/2 + R2Sin(|)2) 


for = 1, 4)^ > tt/2 

Finally, the perimeter of the contact zone for the i 
segment is given by 


(A-60) 


th 


il^(i) = 2 dx 


for = 0 


(A-61) 


and 


£^(i) = B [2 (d/2 + R2Sin(j)2) 

- R2Sin(j)2 + R^Sin4>^J 


(A-61) 


The values of V. , A , and A for the full trunk and 

t c' c g 

cushion are obtained by subtracting the r values from the i 
values for each segment and summing them over all the segments. 


V, 


4 (N+M) 

[V(l) - v^^(l)] 

i=l 


(A-35) 


4 (N+M) 
i=l 


(A-62) 
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(A-63) 


4 (N+M) 


V 




+ V. 


i=l 


and 


4 (N+M) 


Y Ki ■ -^gr 


i=l 


The values of A^^, and A,.^^ are 

by adding the values for each segment. 


tc 


4 (N+M) 


Y ''tci'i' 


i=l 


ta 


4 (N+M) 

E 

i=l 




cn 


4 (N+M) 

E 

i=l 


A . (i) + A (i) 
cnx cnr 


tac 


4 (N+M) 




i=l 


-tcc 


4 (N+M) 




1 = 1 


(A-64) 


obtained 


(A-65) 


(A-66) 


(A-67) 


(A-68) 


(A-69) 
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Center of Pressure 


The distance of the centers of pressure of each segment from 
the center of the cushion are required in order to estimate the 
torques acting on the ACLS. The positions of the centers of 
pressure depend on whether or not the segment is in ground con- 
tact as illustrated in Figure 55. 


Without Ground Contact - Since the pressures inside the 
trunk and the cushion are uniform, the pressure centers coincide 
with the respective centroids of the projected area. The center 
of pressure distances X^j^(i) and Z^j^(i) are: 


^s = 1 



>=ch<i> 

= - L^/2 - 

(d/2 + R 2 Sin()) 2 ) ^2 Cos6 (i) 

(A-70) 

where 





^2 

1.333 

6 

Sin(B/2) 





= Z (i) 

cx 2 


(A-71) 

where Z 

cx 

(i) is evaluated from 

Equation (A-14) . 


1=2 

s 

=^ch<i> 



(A-72) 





(A-73) 
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where X (i) and 

wX 

"cx 

(i) 

are obtained 

from 

Equations 

(A-15) 

and 

(A-16) respectively. 







^s = 3 








='ch'i> 

= 

^cx 

(i) 




(A-74) 

"ch<“ 

= 

"cx 

(i)/2 




(A-75) 

where X (i) and 
cx 

z 

cx 

(i) 

are obtained 

from 

Equations 

(A-17) 

and 


(A-18) respectively. 



X^j^(i) = Lg/2 + (d/2 + R 2 Sin<j> 2 ) Sin6 (i) (A-76) 

^oh<i> = ^ox<i> ®2 

where Z^^(i) is obtained from Equation (A-20) . 



Xch(i) = Lg/2 + (d/2 + R2Sin(f.2)Cos6 (i) (A-78) 

"oh<i> = »2 

where Z (i) is evaluated from Equation (A-22) . 
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6 


II 

•H 

O 

X 

X (i) 
cx 


(A-80) 


Z (i)/2 
cx ' 


(A-81) 

where X (i) and Z 

(i) are obtained 

from Equations (A-23) 

and 

(A-24) respectively. 

1=7 




s 


X (i) 
cx ' 


(A-82) 

"oh<i) = 

Z (i)/2 
cx ' 


(A-83) 

where X (i) and Z 

cx cx 

(i) are obtained 

from Equations (A-25) 

and 

(A-26) respectively. 

1=8 

s 

===h<i> = 

- [Lg/2 + (d/2 

■f R2Sincj>2) Sin6(i)B2j 

(A-84) 


»2 


(A-85) 


where Z (i) is obtained from Equation (A-28) . 


with Ground Contact - The center of cushion pressure dis- 


tances X^j^(i) , Z^j^(i) and the center of trunk pressure distances 
Z^j^(i) for the i^^ segment in ground contact are 


Xt]^(i) , 
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1 


= - L^/2 - 32 (d/2 + R2Sin<J)2 

- R 2 Siri(ti 2 ) Cos6 (i) (A-86) 

X^j^(i) = - Lg/2 - X ^2 Cos6(i) (A-87) 

where 

_ 4 Sin(B/2) 

^x2 3 3 

R^ = d/2 + R2Siri(i)2 + R^SincJi^ (ij)^ _< tt/2) 

R^^ = d/2 + R2Sin(ti2 “ R2Sin<j)2 






32 (d/2 + R2Sin<j)2 - R2Sin(#i2) 

Sin6 (i) 

(A-88) 



= 

X ^2 Sin5(i) 


(A-89) 

= 

2,3 





s 








= 

X (i) 
cx 


(A-90) 


’'tk'i' 

= 

X (i) 
cx 


(A-91) 

where 


obtained from Equation (A-15) for 

"s = 2, 

and 

from 

Equation (A 

-17) 

for Ig = 3. 



■ 


= 

(d/2 + R2Sin4>2 ~ R2Sin(|)2)/2 


(A-92) 
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(A-93) 




= d/2 + R2Sin<|)2 


6,7 


+ (R^SiiKt)^ - 'R^SLn^ /2 




= Lg/2 + ^2 + R2®^”‘t’2 


- R 2 Sin(j) 2 ) Sin6 (i) 


(A-94) 




L^/2 + X ^2 Sin« (i) 


(A-95) 


2ch<i> 


= 32('^/2 + ~ Cos6(i) (A-96) 




X^2 Cos6 (i) 


(A-97) 




Lg/2 + ^2 (d/2 + R2Sin(|>2 


- R 2 Sin(|) 2 ) Cos6(i) 


(A-98) 


X^k(i) 


Lg/2 + X ^2 Cos6(i) 


(A-99) 




+ R2Sin4>2 - Sin6 (i) (A-lOO) 




~^x2 


(A-101) 


='oh<i> 


X (i) 
cx ' ' 


(A-102) 
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= X^^(i) (A-103) 

where X (i) is given by Equation (A-23) for 1=6, and by 
Equation (A-25) for Ig = 7. 

Z^j^(i) = - (d/2 + R2Sin<j>2 - R2Sinc|)2)/2 (A-104) 

Z^j^(i) = -d/2 - R2Sin(|)2 - (Rj^Sincj)^ 

- R2Sin(j)2)/2 (A-105) 




X 

Z 


tk 

ch 


(i) 

(i) 


z 


tk 


(i) 


It should be 
equations is tt/2. 


= - Lg/2 - 02 

- R 2 Sin<t) 2 ) Sin6(i) 

= - Lg/2 - X ^2 Sin5(i) 

= - 02 

- R 2 Sin<{)^) Cos6(i) 

= - X^2 

noted that the upper bound on ^ 


(A-106) 

(A-107) 


(A-108) 

(A-109) 


in these 
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The Flow Model 


The flow through system is shown in Figure 56. The 
through the upstream fan orifice is given by 


Qf ^af ^af 


-2P 


af 


where 


= volume flow through the fan 


= orifice area, atmosphere to fan inlet 


= fan inlet orifice discharge coefficient 


’af ~ inlet pressure (negative, gage) 


p = air density 


The fan pressure rise, P^, is given by 


P^ = P - P 4= 
f p af 


where 


Pp = plenum pressure (gage) 


flow 


(A-110) 


(A-111) 
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The fan flow is obtained from the fan characteristic, 
is a function of the fan pressure rise. 


Qf = f(Pf) 


The remaining flows are found as follows. 


where 


2 P 


Q = A C / — 
'^pa pa pa \ p 


pa 


plenum-atmosphere flow rate (excluding relief 
valve flows) 


= plenum to atmosphere orifice area 


C = discharge coefficient for A 
pa ^ pa 


where 


where 


Qr: = Qx: - Q 

fp ^f ^pa 


= fan-to-plenum flow rate 


= Q + Q ^ 

^fp '^pc ^pt 


Q = plenum-to-cushion flow rate 
pc 


Q . = plenum-to-trunk flow rate 

pt 


and 


(A-112) 


(A-113) 


(A-114) 


(A-115) 


167 



2 CP - P ) 

Q = A ^ C J 2 — 

^pc pc pc \ p 


(A-116) 


where 


= cushion pressure (gage) 


A = plenum-to-cushion orifice area 
pc ^ 


C = discharge coefficient for A 
pc ^ pc 


(Q. 


V 


0, since in static operation, the pressure relief 
valve is closed) 


Q ^ = A ^ C . 

^pt pt pt 


2(P - P ) 

P t 


(A-117) 


where 


P^ = trunk pressure (gage) 


A = plenum-to-trunk orifice area 
pt 


C ^ = discharge coefficient for A . 

pt ^ pt 


^pt ^tc ^ ^ta 


(A-118) 
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where 


Q. = trunk- to-cushion flow rate 

tc 

= trunk-to-atmosphere flow rate 



where 


(A-119) 


tc 


tcc 


trunk-to-cushion orifice area 


orifice area on cushion side of the trunk-ground 
contact 


C. = discharge coefficient of the trunk holes on the 
cushion side 


The factor 2/3 comes in Equation (A-119) due to the trian- 
gular pressure profile assumed in the trunk-ground contact zone 
(see Figure 17) . Since the pressure on the cushion side of the 
contact zone is linearly increasing from P to P . , the flow is 2/3 
of what it would be if the pressure was unform and equal to P^ • 

A^^^ is zero if there is no trunk-ground contact. 


Q 


ta 



(A-120) 


where 


A^^ = trunk-atmosphere orifice area 
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A = orifice area on atmosphere side of trunk contact 

C. _ = discharge coefficient of the trunk holes on the 
atmosphere side 

Just as in Equation (A-119) , the factor 2/3 comes in 
Equation (A-120) due to the pressure profile assumed in the 
trunk-ground contact area. Since the pressure on the atmosphere 
side of the trunk contact zone is linearly reducing from to 0 
(gage) , the flow is two-thirds of what it would be if the pres- 
sure was uniform and equal to zero. 


is zero if there is no trunk contact with the ground. 


Q = Q + 
ca '^pc tc 


(A-121) 


where 


= cushion-to-atmosphere flow rate 


Q 


ca 



(A-122) 


where 



cushion gap area 


C = coefficient discharge for A 

g g 
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The Force Model 


The vertical force developed by the cushion is given by 

F = (P A + P.A ) Cose C 0 S(j) (A-123) 

cn c c t cn e ^e 

where 

P^ = cushion pressure, gage 

A = cushion area 
c 

P^ = trunk pressure, gage 

A = trunk-ground contact area 
cn ^ 

= Eulerian roll angle 
<f)g = Eulerian pitch angle 

The torque developed about the cushion center by the cushion 
pressure force and trunk contact force about the vehicle x axis, 

T , and vehicle z axis, T , are: 

XI z 

4 (M+N) t 
i=l ( 
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4 (M+N) ( 

T = 1 P Fa . (i) - A (i)l [x 

nz / j I c L Cl cr J L ch J 

i=l ( 

+ F-^ • (i) + A (i)l [x., (i)J 

t L cm cnr -* >- tk 

Under equilibrium conditions 


F 

= M^g 

(A-126) 

cn 

a 

T 

nx 

= -F^M^g/ (Cos(j)^ Cos0^) 

(A-127) 

T 

nz 

= C^M^g/(Cos<j)^ Cos0g) 

(A-128) 


(A-125) 


Under the equilibrium loading, the aircraft orients itself 

at a particular For a given value of X^^, the 

variables Y , i> and 0 uniquely define the aircraft and ACLS 
eg e e j 

position (Z and ijj are zero) . The area and distance variables 
^ eg e 

needed to evaluate Equations (A-123) , (A-124) and (A-125) are 
found from the relationships derived in the previous sections, 
and can be expressed as 


A . (i) 

Cl 

II 

!-■ 

^^cg' 



(A-129) 

A (i) 
cr 

= ^2 




(A-130) 



'%■ 


^e) 

(sum of A . (i) and (A-131) 


= ^4 

•’'eg' 



(A-132) 


= ^5 



^e) 

(A-133) 
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(A-134) 


and 


oh<i> 

= ^6 (^cg' 

6 , (j> ) 

e ^e 


= ^7 (^cg' 

e , (b ) 
e' ^e 


(A-135) 


Also the orifice areas required fay Equations (A-119) , 
(A-120) and (A-122) depend on Y / 6 , (j) and can fae expressed 
as 


\c 

= ^8 


■e- 

(D 

CD 

(A-136) 

^ta 

= Fg 


0,6) 
e ' ^e 

(A-137) 

A 

g 

^10 


0 , (j> ) 

e ^e 

(A-138) 


For the static solution, it is necessary to solve Equations 
(■^“110) to (A-138) simultaneously to determine the equilibrium 
aircraft position and ACLS pressures and flows. 
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THE DYNAMIC MODEL 


The dynamic behavior of the ACLS is determined from the 
simultaneous solution of the state equations describing the body 
dynamics and fluid mechanics of the system. 


The Force Model 

The following assumptions are made in deriving the equations 
of motion. 


a) The aircraft does not experience any rotational 

motion about the vehicle y axis, i.e., there is no 
yaw motion in the vehicle coordinate frame 

(ip = ip = ip = 0) . 


b) Forces acting along the vehicle z axis are zero 

(F^ = 0). 


C) 


The aircraft CG stays in the inertial xy plane. 


1 . e , 


ze 


= 0 and Z 


eg 


= Z 


eg 


= Z 


eg 


0 . 


d) The component of the vehicle x axis force along the 

inertial vertical axis is negligible compared to the 

corresponding component of the vehicle y axis force, 

i.e., F Sin<f>„ Cos0 << F Cos(j) Cos0 . 

' X ^e e y e e 


e) The forward motion of the aircraft CG (i.e., the 

motion along the inertial X axis) occurs at a constant 
deceleration which is determined from the braking 
coefficient of the ACLS. 


f) All forces act along the appropriate vehicle axes and 
not along the inertial axes. This comes about because 
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I 


the ACLS model has been set up such that the ground 
underneath any particular trunk segment is considered 
parallel to the hard surface. 

Forces and Torques 

Vehicle y-axis Forces - The dynamic forces along the vehicle 
y axis consist of 


a) The cushion pressure force, F 


cp 


F = A 

cp c c 


(A-139) 


b) The trunk contact force during trunk-ground contact 
tp 


^tp ^t ^cn 


(A-140) 


c) The aerodynamic drag force, 


'df = - Apj, p V /|V| 


(A-141) 


where V, the velocity component along the vehicle 
y axis is 


V = X (Cos(|)„ Sin0„ Sinij^ - CosJ; Sind) ) 
eg e e ’^e e e 


^cg (Cos<f)^ Cose^) 


d) The trunk damping force during trunk-ground contact, 

F 4. 

ct 
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ct 


4 (M+N) 

^ segment is in 

i=l ground contact 

(A-142) 


0 if the segment is not in ground contact. 

v/here the velocity of the trunk segment with respect to ground, 

V^(i) is given by the addition of the component of along the 

vehicle y axis and the component along the vehicle y axis of the 

relative velocity of the trunk segment center with respect to the 

CG, V (i) . The first component is Y Cosd) Cose . The second 
r eg e e 

component is 


V (i) = 0 ) X r (i) 

r ' ' I y axis 


where 


->■ 

u 


(()U + 0u 

^ Z X 


?(i) = [x^j,(D - cj H- [-Gg - 5^ 




and u , u and u are unit vectors along the vehicle x, y and 
X y z 

z axes. After substitution 


' [^cg + * Kk<^> ' ""c J 


- ® [^k<"> - 


(A-143) 
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The total force along the vehicle y axis is 


F = F + F^ + F,- + F . 
cy cp tp df ct 


(A-144) 


and total force along the inertial Y axis is 


F = F Cos(|) Cos0 
y . cy ^e e 


(A-145) 


Vehicle x-axis Torques - The torques acting about the air- 
craft CG on the vehicle x axis consist of : 


a) The cushion pressure torque, 

T = - P fz , (i) - F^1 Ta . (i) - A (i)1 (A-146) 

cpx c L ch fj L cx cr -I 

b) The torque due to trunk contact, 

"^tpx " “ ^t “ ^fl [^cni^^^ "*■ ^cnr^^^^ (A-147) 


c) The torque due to aerodynamic drag, 


'dfx 


= -F C 

df ^fz 


(A-148) 


where Cj^ is the distance along the vehicle z axis of the center 
f z 

of aerodynamic drag from the CG. is positive if the CG is 

left of the drag center. 


d) The torque due to trunk damping, T 


4 (M+N) 


■ tx 


= - E [ 




] "'ct 


(i) 


i=l 


(A-149) 
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The total torque along the vehicle x axis, T^, is then given 


by 


m 5=rp 4.T +T +T 

■^X "^cpx "^tpx "^dfx ■‘^tx 


(A-150) 


Vehicle z-axis Torques - The torques acting about the 
aircraft CG about the vehicle z axis consist of: 


a) The cushion pressure torque, T 


cpz 


’cpz 




b) The torque due to trunk contact, T 


tpz 


■ tpz 


'P^. rx.v(i) “ C "1 Fa . (i) + a (i)1 (A-152) 

t L tk cj L cm cnr J 


c) The torque due to aerodynamic drag, T 


df z 


T = F C 

^dfz -^df ^fx 


(A-153) 


where is the distance along the vehicle x axis of the center 

of drag from CG. is positive if the CG is behind the 

drag center. 


d) The torque due to trunk damping, T 


tz 


4 (M+N) 


tz 


IZ - ''c] 


(A-154) 


i=l 
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e) The torque due to ground friction, T 


fz 


4(M+N) 

- [a ■ (i) + A (i)"l 

/ - t L cm cnr J 

i=l 


■f z 


Ty , (i) + G , , - 
L gh' ' gj ’ eg 


1 ; X ^ 0 

tJ eg ^ 


(A-155) 


; X =0 
eg 


The total torque along the vehicle z axis, T^, is then given 


by 


T = T + T. + T,j; + T^ + Tj; 
z cpz tpz dfz tz fz 


(A-156) 


The Equations of Motion 


The equations of motion give the relation between the forces 
and torques acting on the aircraft and its linear and rotational 
motion. 


Linear Motion - The equations of linear motion are derived 
by applying Newton's second law in the inertial coordinate frame. 
The vertical motion of the aircraft CG is given by: 


M Y 
a eg 


F 

y 


-Mg 


(A-157) 


where F^ is the inertial vertical component of the vehicle y axis 
force. The equation of motion along the inertial x axis (forward 
motion) is simply 
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-y 


(A-158) 


where y is the specified deceleration rate. 

Rotational Motion - If the total torque acting about CG of 
the aircraft is denoted by the vector T, then 

T = dH/dt (A-159) 

where H is the angular momentiam vector about the CG in the 
inertial frame. 

The rate of change of momentum, dH/dt, can be expressed as 
the sum of two components; one describing its rate of change 
relative to the vehicle frame and the other accounting for the 
rotation of the vehicle frame relative to the inertial frame. 

Thus 


dt 




f. 


(A-160) 


where w is the angular velocity vector of the vehicle frame 
(relative to the inertial frame) and is given by 


= 0U + ljjU + d)U 

X y z 


(A-161) 


where u^. 


u^ and u^ are unit vectors along the x, y and z axes 
of the vehicle frame. 


The angular momentum vector H and torque vector T can also 
be divided into instantaneous components along the vehicle axes 
as follows 
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(A-162) 

(A-163) 

(A-164) 
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Expressing Equations (A-159) and (A-160) in matrix form. 




T 


z- 


I 

I 

X xy 

xz 

I 

I 

xy y 

yz 

I 

I 

zx yz 

z 


0 

. i|j 

- <t> - 


- 0- 


ri 

I 

I 1 



X 

xy 

xz 


X 

I 

I 

I 



xy 

y 

yz 

6 


,I 

I 

I . 

w r M 

zx 

yz 

z 



(A-165) 


Equation (A-165) gives three simultaneous equations which 

• • • • • • 

can be solved to obtain the rotational accelerations 9 , tp and 

in terms of the angular velocities and torques. Substituting 

• • • • • • • 

^ = ^ = ^ = 0 (no yaw motion) and solving for 0 and ((> gives 

II - I^ 

Z X zx 


<p = 




I 0+1 
xy yz 


♦>] - "zx 


(A-167) 


II - I 
Z X zx 


Angular Coordinate Transformation 


The angular position of the aircraft in the inertial frame is 
defined by the Eulerian angles 0^, ip^, and (j)^. To obtain a 
specified angular orientation, the aircraft must be rotated in 
the following sequence. First, rotation through angle \p^ about 
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the inertial Y axis. Second, rotation through angle 0^ about 
an intermediate axis (i axis) , which coincides at this point only 
with the vehicle roll axis (x axis) . Finally, rotation through 
angle about the vehicle z axis. With this definition, the 
angular velocity vector can be expressed as 


CO 


♦e “y + 


0 u. + 6 u 
ex e z 


(A-168) 


By solving for the projections of the vehicle axes on the 
Euler angle rotation axes the following relationships can be 
obtained . 



Uy 

-> 

^3 

Cos0 

e 

u Sin0 
z e 

(A-169) 

where 

->■ 

= u 

y 

COScb 

^e 

u Sind) 

X ^e 

(A-170) 

and 

^i 

->■ 

= u 

X 

Cosd) 

- u SiniJ) 
y ^e 

(A-171) 


From Equations (A-161) and (A-168) to (A-171) the relation 
between vehicle frame and Euler angle velocity components is 
found as follows 


Sind) Cos0 + 0 CosiJ) 
^e e e ^e 

(A-172) 

Cosd) Cos0i - 0 Sin4) 
^e e e e 

CA-173) 

• 

- i|) Sin0 
^e e 

(A-174) 
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Substituting t}i = 0 (no yaw) the above equations can be 
solved to give 

• • 

6 =0 Coscj)^ (A-175) 

© © 

• • • 

(}> = (f) + 0Sin<j) Tan0 (A-176) 

• • 

= 0Sin<j) Sec0 (A-177) 

e e e 

The above three equations are the required transformation 

• • 

equations which convert velocities in the vehicle frame (0, (j> 

• • • • 

and ip) to those in the Eulerian frame (9^/ <P^ and '1'^) • The 
Eulerian angular velocities can then be integrated to uniquely 
fix the instantaneous orientation of the aircraft. 


The Flow Model 


The flow model determines the variations in pressures and 
flows as a function of time. There are three parts to the flow 
model: the fluid chambers (i.e., plenum, cushion and trunk), 

the pressure relief valve, and the fan. The principal assump- 
tions of the flow model are as follows. 

a) The flow through all orifices is one-dimensional 
and quasi-static, i.e., the pressure in the plane of 
the orifice is uniform, and the unsteady state terms 
in Bernoulli's equation are small compared to the 
change in velocity head. 

b) The flow through the orifices is incompressible , i . e . , 
the pressure drop is small compared to the total 
pressure, and the air density is constant. 
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c) The pressure and volume changes of the air during 

expansion and compression in the various fluid 

chambers are governed by a polytropic relationship, 
k 

i.e., pv = const.. 


The Fluid Chambers 


Plenum - From the polytropic pressure-density relation 




= constant 


(A-178) 


Taking time derivatives. 


dP 

dt 


k(P + P ) 

p a . ^ 

p dt 


(A-179) 


Conservation of mass in the plenum requires that 


^ < 0 ''■p> 


= P - P Q 


pa 


-pQ -pQ.-pQ 
pc pt V 


vzhere = dynamic fan flow 

From Equations (A-179) and (A-180) , 
d 


dt ^p 


k(P + P ) 
P 




(Qfj^ - Qp^ 


(A-180) 


- Q - Q a. - Q ) 

'^pc ^pt V 


(A-181) 
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Cushion - The continuity equation for the cushion is similar 
to that of the plenum but with an additional term to include the 
rate of change of cushion volume due to motion. 


dt 


k(Pc + P^) 
V_ 


Q + - Q 

pc '^tc ca 


dV 

c 

dt 


(A-182) 


where 



rate of change of cushion volume. 


Trunk - The continuity equation for the trunk is similar to 
that of the cushion 



k(Pt + p^) 


^pt ^tc ~ ^ta 


dt 


(A-183) 


where 



rate of change of trunk volume. 


The Pressure Relief Valve 


A schematic diagram of the pressure relief valve is shown 
in Figure 57. The valve is modelled as a second order mass- 
spring-dashpot system. The equation of unrestrained valve 
motion (i.e., the valve motion if stops were absent) is found 
from the relationship below. 


m X + Zx + kx 

V V V V V V 


(P - P , ) A. 
P pb f 


(A-184) 
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Valve 

Stem 


Atmosphere 



Figure 57. Schematic Diagram of Pressure Relief Valve 
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The above equation neglects the aircraft acceleration in 
comparison to the valve stem acceleration. 

Due to the presence of stops , there are two constraints on 
valve motion. 


If 

X 

£ 0 , and 

P < P w 


V 

• • 

P Pb 

Then 

^v 

= X = 

V 

0 

If 

X 

> X and 

P A^ > P Aj; + k X 


V 

• 

— a 

p f pb f V a 

Then 

X 

= X = 

0 


V V 


where is the range of valve motion between stops . The above 
constraints ensure that the valve motion does not exceed the 
stop limits. 

The vent area for a given valve displacement x^ is found 

from 


A 


V 



(A-185) 


where is an equivalent valve periphery which can be a constant 
or be a function of the valve displacement depending on the valve 
design. The flow through the valve, which depends on the vent 
area A^, is given by 


Q 


V 



(A-186) 
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where C is the discharge coefficient of the valve and is 
the number of (identical) relief valves. 


The Fan 


The dynamic equations for the fan are 


dQ 


fx 


dt 


+ P - P 
f ar £ 


(A-187) 


where is the inertance of the fan and is the dynamic fan 

f low , and 


^f (A-188) • 

where the function represents the static fan characteristic for 

positive and negative flow rates. The fan volume, represented 

by the fan capacitance is lumped together with the plenum volume, 

so that V in Equation (A-181) includes the added volume due to 
P 

the fan. 
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